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SUMMARY 


Newly emerging aerospace technology points to the feasibility of sustained hypersonic flight De- 
signing a propulsion system capable of generating the necessary thrust is now the major obstacle. First- 
generation vehicles will be driven by air-breathing scramjet (supersonic combustion ramjet) engines. Be- 
cause of engine size limitations, the exhaust gas leaving the nozzle will be highly underexpanded. Con- 
sequently, a significant amount of thrust and lift can be extracted by allowing die exhaust gases to ex- 
pand along the underbody of the vehicle. Predicting how these forces influence overall vehicle thrust, lift, 
and moment is essential to a successful design. This work represents an important first step toward that 
objective. The UWIN code, an upwind, implicit Navier-Stokes computer program, has been applied to 
hypersonic exhaust plume/afterbody flow fields. The capability to solve entire vehicle geometries at hy- 
personic speeds, including an interacting exhaust plume, has been demonstrated for the first time to the 
author’s knowledge. Comparison of the numerical results with available experimental data shows good 
agreement in all cases investigated. For moderately underexpanded jets, afterbody forces were found to 
vary linearly with the nozzle exit pressure, and increasing the exit pressure produced additional nose-down 
pitching moment. Coupling a species continuity equation to the UWIN code enabled calculations indicat- 
ing that exhaust gases with low isentropic exponents ( 7 ) contribute larger afterbody forces than high -7 
exhaust gases. Moderately underexpanded jets, which remain attached to unswept afterbodies, underwent 
streamwise separation on upswept afterbodies. Highly underexpanded jets produced altogether different 
flow patterns, however. The highly underexpanded jet creates a strong plume shock, and the interaction of 
this shock with the afterbody was found to produce complicated patterns of crossflow separation. Finally, 
the effect of thrust vectoring on vehicle balance has been shown to alter dramatically the vehicle pitching 
moment. 
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PREFACE 


Hypersonic flight represents the newest frontier in man’s quest to fly higher and faster. Unique among 
aerospace vehicles, hypersonic craft will roam the diffuse boundary between conventional atmospheric 
flight and space flight; they will bridge a gap in the increasingly traveled path to space. But to do so 
will require quantum technology leaps in several aerospace disciplines. Scramjet propulsion, at the fore 
of hypersonic research activity, provides the impetus for computational aerothermodynamicists to address 
this problem. Applying their approach to develop a predictive method for hypersonic vehicle flow fields 
will facilitate the design process at all stages. Progress toward this goal is being pursued along several 
fronts: geometry definition and grid generation, flow modeling, algorithm development, code calibration, 
detailed component studies, and configuration studies. 

Within the context of scramjet exhaust plume flows, this work attempts to contribute to each of these 
categories, albeit in varying degrees. By encompassing some specific and general aspects of plume flows, 
the investigation offers insight into parameters that govern the problem. It also reflects the current status 
of predictive capabilities for these flows. The following pages describe this investigation in its entirety. 
Relatively new aspects of the problem to computational methods, such as coupling a species continuity 
equation or implementation of the total-variation-diminishing (TVD) scheme, are presented in some detail 
and are referenced for further reading. However, other aspects standard to computational fluid dynamics 
are addressed only insofar as they impact the present discussion: the Baldwin-Lomax turbulence model, 
for example. 
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CHAPTER ONE 


INTRODUCTION 
1.1 Background 


Hypersonic flight is a vital link in the advancement of transportation within and beyond the earth’s 
atmosphere. Whether it requires the relatively short exposure of a reentry vehicle, or the prolonged cruise of 
a hypersonic transport, surviving the aerothermodynamic environment characteristic of this flight regime 
constrains the design of a wide class of aerospace vehicles. This investigation studies the details of an 
important determinant of one type of hypersonic vehicle ’s performance: impingement of the engine exhaust 
plume on the underbody. In order to appreciate the parameters governing the design of a hypersonic vehicle, 
it is helpful first to classify the hypersonic flight condition. 

The boundary between “supersonic” and “hypersonic” flight is somewhat vague. Certainly not less 
than Mach 3, nor greater than Mach 10, the dividing line defies straightforward determination of the point at 
which flight conditions are first considered hypersonic. In fact, the distinction is more correctly associated 
with the energy of the approaching air (with respect to the reference frame moving with the flight vehicle). 
Because the air is brought to a rest at the surface of the vehicle, the stagnation energy provides a meaningful 
measure of the thermal environment surrounding the aircraft. For an ideal gas, the stagnation energy has 
a quadratic dependence on the flight Mach number. Thus, at subsonic speeds, the stagnation energy is a 
small percentage above the static energy, while at supersonic speeds, the stagnation energy rises rapidly. 
In the Mach 3 to Mach 10 interval, the stagnation energy increases from a level about double the static 
energy to nearly ten times that level. The emergence in this interval of flow energy as a prominent feature 
sets hypersonic flow apart from supersonic. 

Whereas pressure loads drive the design and analysis of supersonic vehicles, high-temperature phe- 
nomena assume central importance for hypersonic flight. As the flight Mach number increases, strong 
shocks can energize the air sufficiently to cause excitation of new internal energy modes and enable chem- 
ical reactions among the component species of air. These phenomena, inconsistent with ideal gas assump- 
tions, significantly alter the shock structure, temperature distribution, and transport properties of air from 
the results of an ideal gas analysis. Hypersonic aerodynamics, then, is not simply an extension of super- 
sonic aerodynamic theory; rather, it is a physical science unto its own. Development of a hypersonic flight 
vehicle requires new technology in materials, thermal protection, aerodynamic design and propulsion to 
address the constraints attendant to the hypersonic flight regime [lj. 

1.2 Motivation 


Research stimulated by federally- sponsored programs such as the Space Shuttle, the National Aero- 
space Plane (NASP), the Aero-assisted Transfer Vehicle (AOTV) and the Aero-assisted Flight Experiment 
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(AFE) has given rise to technological advancements in several critical areas. Experts now suggest that sus- 
tained hypersonic flight is not only possible, but also practical [2,3], The promise of reducing trans-Pacific 
flights from a grueling 16 hours to about 2 hours in a Mach 14 airliner has piqued commercial interest. 
Furthermore, airbreathing scramjet engines, which benefit payload capacity dramatically by eliminating 
the need to carry oxidizer, make numerous aerospace missions feasible. However, before any of these 
goals are realized, many technological questions remain to be answered. 

1.3 Problem Definition 


Propulsion is the primary technical problem in designing an air-breathing hypersonic vehicle [4], At a 
cruise speed in excess of Mach 6, few practical power plant choices exist. The hydrogen-powered scramjet 
engine will likely be employed for the first generation of hypersonic vehicles. As suggested by the name, 
the scramjet (supersonic-combustion ramjet) maintains supersonic flow through the combustion chamber. 
The hydrogen fuel mixes and bums rapidly enough with atmospheric oxygen to contain the process within 
a reasonably-sized engine. However, scramjets can operate only at flight Mach numbers above about 4, so 
alternate propulsion systems are required to accelerate the vehicle from takeoff to cruise speed. Undoubt- 
edly, the hypersonic cruise vehicle will require multiple propulsion systems, for example, an afterburning 
turbojet for acceleration from rest to low supersonic speed, a ramjet for supersonic acceleration, and finally 
a scramjet for cruise [5]. 

Aside from the obvious logistical problems of installing and operating several engines on the same 
vehicle, researchers have yet to prove conclusively that an installed scramjet can produce net thrust For 
an airbreathing engine, net thrust represents the difference between gross nozzle thrust and inlet ram drag. 
Preliminary studies indicate that a Mach 12 hypersonic cruise vehicle weighing 500 , 000 lb could experi- 
ence ram drag as high as 1.95 million pounds, requiring up to 2.1 million pounds of gross thrust [6], The 
thrust margin, or percentage of gross thrust in excess of ram drag, in this instance is only 7.7%. Needless 
to say, small deviations from design conditions could easily lead to inadequate engine performance. 

Maximizing thrust and minimizing drag necessitates using the entire underbody to process the mass 
flow through the engine. Such an “airframe-integrated” propulsion system takes advantage of the precom- 
pressed forebody flow upstream of the engine inlet. Ideally, the inlet captures only air that has traversed 
the vehicle’s bow shock. To design an inlet that captured air outside the bow shock would lead not only to 
unacceptably nonuniform inlet entrance conditions, but also to a severe inlet drag penalty attributed to the 
inlet “bow” shock. Upon entering the scramjet inlet, the precompressed capture air is further compressed 
in a series of ramps before reaching the combustion chamber. The combustion chamber consists simply of 
a diverging duct with fuel injectors installed in the sidewalls and the inlet struts. Hydrogen fuel possesses 
the advantage of a high flame speed, which reduces the length of the combustion chamber. However, it 
is also a comparatively low density fuel and hence constrains the vehicle size to accomodate a large fuel 
volume. How completely the fuel will be burned within the combustion chamber is unresolved. 

Airframe integration will also come into play in expanding the exhaust gases. The high energy gas 
exiting the nozzle will be significantly higher than ambient pressure, but installing a nozzle to expand this 
gas completely is impractical. Not only would it be extremely heavy, but it would also create unnecessary 
drag due to the additional wetted area. In fact, very little of the exhaust gas expansion will be carried out 
internally. Instead, the exhaust gas will expand externally, washing the afterbody with an underexpanded jet 
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flow. In this region, the flow will bear similarities to both a diverging nozzle flow and an underexpanded 
free jet. Careful contouring of the afterbody shape is critical to extract all the available thrust from the 
engine. Equally important, the high pressure gas impinging on the underbody will contribute to the lift and 
pitching moment of the vehicle, and these effects must be accounted for in locating the engine, sizing the 
wing planform, and assessing control surface requirements. 

Clearly, design of the hypersonic vehicle’s propulsion system will play a significant role in the over- 
all design of the vehicle itself. A successful design requires detailed understanding of the flow physics 
throughout the propulsion system, along with the capability to predict how propulsion system performance 
is influenced by the geometry of the vehicle. A brief review of prior and ongoing work will help place the 
present work in perspective. 


1.4 Literature Review 


Integration of airframe and propulsion systems in transonic and supersonic flight has been an area 
of active research over the past decade. Investigations of inlet flows using Navier-Stokes computational 
fluid dynamics (CFD) codes have been performed for wing/nacelle and fuselage/inlet geometries, and 
mixed compression ramp inlets [7-9]. For a review of current literature on this subject, see reference 
[7]. A proliferation of work has also come forth recently concerning rocket and missile base flows at 
supersonic speeds. These include axisymmetric Euler methods [10], axisymmetric Navier-Stokes zonal 
methods [11], axisymmetric thin-layer and full Navier-Stokes methods [12], and three-dimensional Navier- 
Stokes methods [13]. However, in comparison to the deliberate progress seen in transonics and supersonics, 
research in hypersonics has experienced a more spasmodic history. 

Spurred by the Apollo program and an ambitious planetary exploration agenda, the late 1960s and 
early 1970s saw rapid progress in hypersonics research. Researchers addressed the problem of scramjet 
propulsion by performing experiments on engine components in high-speed wind tunnels, and by applying 
proven analytical tools such as the method of characteristics and Newtonian flow theory [14-19]. But 
research activity declined sharply in the early 1970s, following the successful completion of the Apollo 
program and preliminary design of the Space Shuttle. Hypersonics research underwent a decade-long 
hiatus, brought to an end within the last few years by the technological requirements of a new generation 
of hypersonic flight vehicles. 

In an effort to make up for lost time, many existing facilities are being brought out of mothballs, from 
shock tubes and ballistic ranges to hypersonic wind tunnels and arcjets. Because of the importance of this 
research to national security, only a few current experimental results in the hypersonic regime are available 
in open literature. Lockman [20] is presently studying the unclassified hypersonic all-body vehicle in the 
NASA Ames 3.5-foot hypersonic tunnel. Lockheed has tested a generic hypersonic vehicle at transonic and 
supersonic Mach numbers in a power-on simulation (private communication from D. Schuster, Lockheed- 
Georgia Co., Marietta, Georgia). Regardless of the availability of such data, ground-based test facilities 
pose an overriding shortcoming: very few of them are capable of simulating the high altitudes, Mach 
numbers, and Reynolds numbers of actual flight conditions. Computational fluid dynamics is poised to 
supply that capability. 

The past decade has been witness to a meteoric advancement in CFD capabilities. A synergis- 
tic coupling of improvements in computer resources and numerical solution algorithms has resulted in 
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progressively faster, more accurate, less expensive flow solutions [21]. MacCormack’s predictor-corrector 
method, now an industry standard, opened the door to Navier-Stokes calculations [22]. The advent of 
practical implicit methods, such as the Beam- Warming algorithm [23], cut computer run times by an or- 
der of magnitude or more, making complex three-dimensional flow solutions possible. By applying a 
diagonalization procedure to the Beam- Warming algorithm, Pulliam and Chaussee cut computer time by 
another order of magnitude [24]. The generality and efficiency of these algorithms has given rise to CFD 
research on problems with unprecedented geometrical and physical complexity. Numerous investigators 
have demonstrated flow solutions for viscous flow about realistic three-dimensional geometries at transonic 
and supersonic speeds [9,25]. 

The biggest obstacle to extending these codes to solve hypersonic flows has been capturing the 
strong shock waves associated with the forward-facing surfaces. Development of flux-split [26] and total- 
variation-diminishing (TVD) algorithms [27-30] has ameliorated this problem substantially, bringing forth 
a generation of robust codes capable of solving hypersonic flows [31-33]. These codes are now being ap- 
plied to entire hypersonic vehicles and to specific components, mainly those associated with the propulsion 
system [34-38]. For the first time, problems are being attempted which model three-dimensional geome- 
tries and real gas flows, either in chemical equilibrium [39,40] or nonequilibrium [41]. At this time, it 
is appropriate to reach beyond code calibration analyses for these flow fields, and address performance- 
related issues associated with scramjet propulsion. 

1.5 Objective 


The goal of the present work is to demonstrate a capability to solve complex nozzle/afterbody prob- 
lems in hypersonics. It is first sought to calibrate the computational results with available experimental 
data, then to analyze the impingement of an underexpanded jet on the vehicle afterbody. This informa- 
tion will contribute to a predictive method for measuring installed scramjet performance, particularly as it 
relates to afterbody geometry and the interacting flow of external fluid and exhaust gases. 

The following pages present an explanation of the theories, methods, and results associated with 
this investigation. First, the governing equations are presented with a description of their assumptions 
and limitations. The numerical methods employed to solve these equations are described, along with 
boundary conditions and the zonal grid methodology. Next, the method is applied to two test cases. First, 
a simplified nozzle/afterbody geometry is studied to gain an understanding of parametric influences on the 
plume/afterbody interaction. Then a more realistic generic hypersonic vehicle is analyzed to examine the 
effect of the exhaust plume on the aerodynamic forces of an entire body. Experimental data calibrate the 
numerical results wherever possible. Finally, some conclusions are drawn about this application of CFD 
to the analysis of hypersonic propulsion systems. 
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CHAPTER TWO 


GOVERNING EQUATIONS AND 
NUMERICAL METHODS 


2.1 Introduction 


This chapter presents the partial differential equations governing fluid flow, die Navier-Stokes equa- 
tions. The thin-layer approximation is invoked, and the equations are Reynolds-averaged to facilitate mod- 
eling turbulence. The equations are then transformed to generalized curvilinear coordinates to simplify the 
boundary conditions and enable efficient use of the computer. 

Because some of the flowfields being considered are composed of binary gas mixtures, a species conti- 
nuity equation is also included. In its most general form, this equation (only one equation is required to trace 
two species) includes terms related to convection, diffusion, production, and destruction of a species. How- 
ever, neglecting the diffusion, production and destruction terms greatly simplifies solving this equation. A 
loosely-coupled approach is applied to solve the system of equations wherein the Navier-Stokes equations 
are solved simultaneously, then the species equation is solved as a separate, scalar equation. Informa- 
tion about the species concentrations is fed back to the Navier-Stokes equations through an equivalent-7 
approach, which is also described in this section. 

An overview of the Osher TVD scheme implemented to solve the Navier-Stokes equations follows. 
This upwind scheme features enhanced stability over conventional central-differencing schemes, and is 
well suited to hypersonic applications. The algorithm has been implemented within a three-dimensional, 
zonal grid approach that is adapted to discontinuous geometrical variations, as occur at the exit plane of a 
scramjet nozzle. The chapter concludes with a discussion of the turbulence model as applied to the flow 
fields of interest here. 


2.2 Governing Partial Differential Equations 


2.2.1 Navier-Stokes equations- In high Reynolds number flows, viscous effects are confined to a thin 
layer near solid surfaces. Viscous derivatives in the stream and crossflow directions are negligibly small, 
while derivatives normal to solid surfaces dominate near-wall phenomena. Retaining only the viscous 
derivatives in the body-normal direction leads to the thin-layer approximation. Transforming the equations 
from Cartesian coordinates to generalized curvilinear coordinates allows solid surfaces to be mapped to 
constant coordinate surfaces, thereby allowing the thin-layer approximation to be invoked for arbitrarily- 
oriented surfaces. The thin-layer Navier-Stokes equations written in generalized curvilinear coordinates 
are 

d r Q+ d{E+ d n F+ d<;d= Re-'d<S (2.1) 
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The vector of conservative variables Q is composed of density (p), x-, y-, and 2 -momentum (pu, pv, and 
pw, respectively), and total energy (e). The pressure is designated by p, and etc. represent the 
transformation metrics. U , V, and W are the contravariant velocity components, defined as 


U = + £xu + £v v + i* w 
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The viscous flux terms are given by 
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In the foregoing equations, density is nondimensionalized by poo. Cartesian velocity components are 
nondimensionalized by Ooo/>/7- and energy by poo- Temperature, a derived quantity, is nondimensional- 
ized by Too. The Reynolds number is symbolized by Re, the Prandtl number by Pr, and 7 represents the 
ratio of specific heats. The dynamic viscosity is symbolized by p. 

Pressure is related to the conservative flow variables Q by the ideal gas equation of state 


p = (7 - 1 ) 



— p(u 2 + v 2 + vr) 


( 2 . 6 ) 
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The metric terms are defined as 
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is the independent variable transformation Jacobian. 
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2.2.2 Species continuity equation- Gas leaving the scramjet nozzles is highly energized, and consists 
of a mixture of fuel, combustion products, and the air that passes through the engine. Whereas the external 
flow is comprised primarily of molecular nitrogen (N 2 ) and oxygen (O 2 ), the nozzle flow of an airbreathing, 
hydrogen-powered scramjet will also contain water (H 2 O) and various intermediate reaction products. 
Indeed, for some flight conditions the flow will still be reacting as it leaves the nozzles. In any event, the 
exhaust flow will certainly possess markedly different thermodynamic properties than the external flow. 
These must be accounted for in the governing equations to properly simulate a scramjet plume flowfield. 


As a first approximation, it is assumed that the chemical composition of the exhaust flow remains 
fixed after leaving the nozzle. This approximation models “frozen” flow, for which the convection time is 
very short compared to the chemical relaxation time. In other words, the exhaust gas traverses the body so 
fast that its composition effectively remains the same from the nozzle exit to the end of the body. 

It is further assumed that the exhaust gases can be lumped together as a single species, and the external 
flow a second species. Thus, the flow consists of two streams of gases that have constant, but different, 
isentropic exponents: coflowing, dissimilar ideal gases. Accounting for mixing processes downstream of 
the nozzles requires determining the concentration of each species throughout the flow field. This is done 
by solving a species continuity equation. In its most general form, the species continuity equation is given 
by 

d(pa) S(pciU + ji.) . d(p(civ + ji,) , d (paw + >,,) 

~dT + — Tx — + Fy + ai K> <29> 

where c* = pi/p is the mass concentration of species i [42]. Ki is the net rate of production of species i 
per unit volume, and (j is , j if , j it ) is the diffusive mass flux vector defined by 
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D{j is the multicomponent diffusion coefficient, and Df is the thermal diffusion coefficient. The quantity 
d ; is the macroscopic gradient vector, given by 
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In the foregoing equations, the n* are the number densities of the species, and n is the overall number 
density of the mixture. The m* are the species molecular weights, and X* represents the external force 
vector acting on the molecules. 


Inspecting these equations reveals that there are four distinct sources of diffusion: (1) diffusion due 
to a concentration gradient, (2) pressure diffusion, (3) diffusion due to external forces, and (4) thermal 
diffusion. When these equations are Reynolds-averaged, a fifth type of diffusion appears: an “apparent 
diffusion” arising from averaging the unsteady convective terms, analogous to the Reynolds stresses that 
appear in the momentum equations. The role of ( 1 ) and (4) were assessed by Ahtye [43], and later (2) 
and (3) by Howe and Sheaffer [44] for the stagnation region of reentry vehicles. The order of magnitude 
analysis of Howe and Sheaffer indicates that, away from shocks and outside the boundary layer, the dif- 
fusive terms are negligible compared to the convective terms in the species continuity equation. A more 
careful analysis is contained in the appendix, so for now it is sufficient to say that the cases under consider- 
ation here feature high flow velocities and comparatively low pressure and temperature gradients, therefore 
convection dominates the solution for species concentrations. 


Eliminating the diffusion terms dramatically simplifies the species continuity equation. Also, by ex- 
panding the differentials, terms corresponding to the mass continuity equation can be omitted, leaving 


dc dc dc dc 

ai + u ai* v a^ +w a I = 0 


( 2 . 12 ) 

When this equation is transformed to generalized curvilinear coordinates, it is written as 

(2.13) 

where r = t, and U, V, and W are the contravariant velocity components, defined in equation (2.3). 


dc dc dc dc 

— + u— + V— + W— = 0 
dr U di dr, ^ 


2.2.3 Coupling the equation sets- The species continuity equation is coupled to the Navier-Stokes 
equations through 7 , the ratio of specific heats. While this quantity is constant for a thermally and calor- 
ically perfect gas, its value depends on the composition of the gas. Kinetic theory [45] predicts that the 
isentropic exponent of a gas is equal to ( n + 2 ) /n, where n is the number of degrees of freedom of the 
constituent molecules. Gases made up of diatomic molecules, such as air, have an isentropic exponent of 
1.4, corresponding to the five degrees of freedom (three translational modes and two rotational modes). 
Having determined the species concentrations by solving the species continuity equation, the equivalent 
7 , designated 7 , is found from 


Cl Cpi + ( 1 Cl) Cp2 
ClC Vl + (1 - Cl) Cvj 


(2.14) 


which amounts to the ratio of mass-averaged specific heats. The ideal gas equation of state remains un- 
changed except that 7 is replaced by 7 : 


p= (7 - l) - jp(u 2 + v 2 + w 2 )^ (2.15) 

Now the numerical methods used to solve the Navier-Stokes equations and species continuity equation are 
discussed. 
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2.3 Numerical Methods 


This section outlines the numerical methods used to discretize the partial differential equations. The 
present work is an outgrowth of the UWIN code, which was developed by Rai as a time-dependent code 
compatible with space- marching (PNS) codes [32], The UWIN code implements the Osher upwind TVD 
scheme in an iterative, flux-split, implicit algorithm to be described below [30]. Upwind schemes offer 
advantages over conventional central-difference schemes for problems in hypersonics. To capture strong 
shocks in hypersonic flow, central-difference schemes require large amounts of user-specified artificial 
smoothing, which tends to degrade the accuracy of the solution. On the other hand, smoothing appears 
inherently in upwind schemes, enabling them to capture very strong shocks without arbitrary parameters. 
Implementing TVD logic further enhances the scheme by preventing new extrema from appearing in the 
solution. TVD schemes have demonstrated the ability to capture shock waves and slip surfaces in one or 
two grid points without oscillations [46]. The species continuity equation is solved uncoupled from the 
Navier-Stokes equations in a fully implicit, upwind, approximately-factored algorithm. 

2.3.1 Osher upwind scheme- A lengthy digression into the mathematical foundation of TVD in 
general, or the Osher scheme in particular, is not warranted here. Instead, the important features of the 
algorithm are described. The reader is recommended to references [27-30] for a deeper analysis. To begin, 
consider the Navier-Stokes equations in curvilinear coordinates, given in equation (2.1). Neglecting the 
viscous terms for the moment leaves 
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The numerical fluxes on the right hand side of equation (2.20) are evaluated with the Osher scheme, which 
for first-order spatial accuracy is given by 




( 2 . 21 ) 


The integrals in the above equation are evaluated on three phase-space subpaths in the interval (t, i + 1 ). 


Let 
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and X is the eigenvalue corresponding to the particular path. For a description of the second-order TVD 
scheme, see reference [30]. 


The Reimann invariants employed in the Osher scheme to solve equation (2.21) are derived for an ideal 
gas. The phase-space integrations in this equation are not exact for variable -7 flows, and produced small 
oscillations in the solution variables in regions of concentration gradients. To ameliorate this problem, the 
isentropic exponent that appears in the Reimann invariants was Roe-averaged: 


y/pili + y^2r7«+i 
\/pi + y/Pi+l 


(2.26) 


This technique was found to reduce the magnitude of oscillations to less than 1% of the nominal value of 
the variables. The UWIN algorithm results in 5x5 block tridiagonal equations along the coordinate lines, 
which are solved by sweeping through the grid in each coordinate direction. 


2.3.2 Upwind implicit algorithm for species continuity equation- The form of the species continu- 
ity equation solved in this study results in a scalar partial differential equation with constant coefficients. 
This type of equation is easily solved with an efficient and robust fully implicit algorithm. One-sided spatial 
differencing is applied according to the local flow velocity, and the equation is approximately factored so 
that it may be solved with scalar tridiagonal inversions in each of the three coordinate directions. Applying 
this method to equation (2.13) results in the finite difference equations 

[(AtU)I + 8^]c" = c n 
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where U, V , and W are the contravariant velocity components as defined in equation (2.3), and 5“ is an 
upwind difference defined as, for example, 


( Ci- Ci—l 
1 


if C7 > 0 


\ 5-(cm-cm) if ^ = ° 
, „ if 17 < 0 

t. c,+i a 


2.4 Ttirbulence Model Considerations 


(2.28) 


Resolving the minute scales of turbulent flow requires an extremely fine grid, and gives rise to a wide 
range of time scales. The resultant stiff equation set would require a prohibitive amount of computer time 
under current computational capabilities. Therefore, modeling turbulence has become standard procedure 
for three-dimensional CFD problems. The algebraic eddy viscosity turbulence model of Baldwin and 
Lomax [47] is implemented in the UWIN code, and is unmodified for these solutions. This method replaces 
the molecular viscosity p, by an effective viscosity (/i + fir) in the Navier-Stokes equations. Near the 
wall, the turbulent viscosity is calculated based on Prandtl mixing-length theory, and farther from the wall 
the calculation is based on the vorticity distribution. For a complete description of the Baldwin-Lomax 
turbulence model, see reference [47]. 

A few comments arc appropriate regarding the use of this turbulence model for the present applica- 
tions. First, the model has been shown to overpredict the turbulent viscosity in regions of moderate flow 
separation, causing the predicted separations to be smaller than those observed experimentally [48]. Also, 
the turbulent viscosity decays to zero far from the wall in this model. This means that the shear layer 
between the jet and external flow will develop only from the molecular viscosity, although laminar shear 
layers are known to become turbulent over a very short distance (lecture notes from B. J. Cantwell, Stan- 
ford University, Stanford, California). However, accurate prediction of turbulent shear layers requires a 
different turbulence model altogether one that models turbulent viscosity in the absence of a wall. Since 
predicting shear layer development was not central to this work, no attempt was made to model them, and 
in any case, the grid point density away from the wall was insufficient to resolve a shear layer. 

Finally, some of the cases presented in the following chapter involve coflowing streams of perfect 
gases. The Baldwin-Lomax turbulence model works equally well for both gases, but the constants that 
appear in the model were derived for constant pressure boundary layers at transonic speeds [47]. Correc- 
tions to these constants to account for hypersonic speed and gas mixtures are necessary, at a minimum, 
to properly account for the turbulent viscosity of the mixing gases. However, as a first approximation, it 
was assumed that the turbulent viscosity for both of the gases can be modeled with the Baldwin-Lomax 
turbulence model in its unmodified form. 
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CHAPTER THREE 


ANALYSIS OF NOZZLE/AFTERBODY FLOW FIELDS 


3.1 Introduction 


The aerodynamics of scramjet propulsion can be described as three distinct components: 1) forebody 
compression, 2) internal fuel injection and ignition, and 3) afterbody expansion (see figure 1). Compressing 
the flow on the forebody is accomplished when the oncoming air traverses the bow shock of the vehicle. 
This serves to create a dense, high pressure flow at the engine inlet, thereby eliminating the need for the 
internal compressors used in turbofan and turbojet engines. Needless to say, rotating compressors would 
be altogether impractical at hypersonic speeds, so external compression becomes a necessity. Ideally, 
forebody compression maximizes the mass flow at the inlet face and produces a uniform flow into the 
engine. The greater the mass flow into the engine, the more fuel can be burned, producing greater thrust. 
But approaching stoichiometric combustion of the fuel depends heavily on uniform flow conditions in 
the combustor. Scramjet inlet designers are challenged with producing a forebody and inlet shape that 
optimizes both mass flow and flow uniformity. 

The internal segment of the scramjet propulsion cycle consists of further flow compression, followed 
by fuel injection and ignition. Injectors are placed on the sidewalls and struts so as not to block flow 
through the engine needlessly. Here, turbulent mixing processes will be relied upon to distribute the fuel 
efficiendy, enabling rapid combustion. How well the fuel mixes with the ingested air is an important 
parameter governing scramjet performance. But the combustion kinetics of hydrogen and atmospheric 
oxygen in a supersonic stream represents a primary question being investigated by scramjet designers: 
how completely will the fuel be burned within the combustor, and what will the flow properties be at the 
nozzle exit plane? Hydrogen’s high flame speed acts in favor of more complete combustion internally, 
and indeed, that is the main reason it is chosen over higher density fuels. However, at the upper limit of 
the flight envelope, say Mach 25 at 150,000 feet, the fuel will doubtlessly be only partly burned by the 
time it leaves the nozzle. Furthermore, much of the energy released during combustion will be absorbed 
by endothermic dissociation reactions of the combustion product (H 2 O). Recovering the energy sunk into 
dissociation reactions depends heavily on catalyzing recombination reactions in the exhaust flow. Attention 
then shifts to the reacting exhaust gas as it flows over the afterbody. 

The final stage of the scramjet cycle, arguably the most important to its performance, is the flow of the 
underexpanded exhaust gases leaving the nozzle. Consider a scramjet-driven vehicle’s afterbody, enclosed 
by a control volume as shown in figure 2. Neglecting upstream influence of the exhaust gas, the force on 
the body resulting from the exhaust gas is given by 

F = — f pVV • dA — [ pnn • dA + f cry ■ dA (3.1) 

Jsi JSi Js 2 

In the above equation, F is the gross thrust, dA is an area element with outward-facing normal, and cry 
is the fluid stress tensor at the wall. Equation (3.1) reveals that gross thrust for the scramjet derives from 
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two sources: 1) the pressure and momentum flux at the nozzle exit plane itself (Si ), and 2) forces arising 
from fluid stress acting on the surface of the afterbody (S 2 ). The fluid stress at the wall consists of viscous 
shear stress acting parallel to the wall, and pressure acting normal to the wall. The shear stresses, which 
are realized as skin friction, mainly contribute to drag. Pressure, acting normal to the wall, produces lift, 
thrust, or drag, depending on the orientation of the surface. By incorporating an upswept afterbody, the 
high pressure exhaust gas generates both lift and thrust. A nose-down pitching moment also results from 
afterbody lift. Since propulsion system lift, thrust and moment account for a significant fraction of the 
overall forces on the hypersonic vehicle, effectively utilizing the afterbody will play an important role in 
obtaining the requisite propulsion system performance. The afterbody design, even the materials chosen 
to construct the airframe, will determine whether the energy in the fluid (kinetic, thermal, and chemical) is 
extracted to generate additional propulsive forces, or is left behind in the wake. 

This chapter opens the investigation into exhaust plume/afterbody flow fields by analyzing a simpli- 
fied geometry, the nozzle/afterbody model. To begin, the model and experimental procedure are described. 
A discussion of grid generation for the numerical problem follows, along with a description of the zonal 
grid scheme. Boundary conditions for this geometry are outlined next. To assess the computational effort 
involved in obtaining the solutions, the solution methodology and code performance are enumerated. Fi- 
nally, numerical solutions are presented for the nozzle/afterbody flow fields. The sensitivity of afterbody 
forces to various parameters is assessed and explained based on the results. The chapter concludes with a 
grid refinement study for the nozzle/afterbody geometry. 

3.2 The Nozzle/Afterbody Model 


Over the past two decades, computational fluid dynamics has advanced rapidly from a fledgling ana- 
lytical tool to a mature, multidisciplinary science. Applications of CFD to subsonic and transonic problems 
are diverse and well validated. Practical aerodynamic research in these flight regimes has made a place 
for CFD to play an important role in the preliminary design of aircraft components such as wings, control 
surfaces, and wing-nacelle fairings [8]. However, expanding the “envelope” of CFD capabilities to a new 
class of problems necessitates validating and/or calibrating the method against experimental or flight test 
data, or at least a well-proven, previously developed analytical method. Doing so instills confidence that 
the numerical model is rendering the correct representation of a physical problem. Calibration serves to 
modify assumptions in the numerical model that may no longer be valid for a new problem. In any case, 
applying a CFD method to one or more representative experiments enables the researcher to: 1) establish 
the capability of the code to solve a new class of problems, and 2) expand the information set associated 
with the experiment, thereby bringing forth a deeper understanding of the results. 

Because this investigation does in fact address flight conditions heretofore unattempted by CFD codes, 
an important first step was to verify the accuracy of the numerical results by comparison with available 
experimental data. To this end, two experiments were identified as appropriate test cases. The first, a 
simplified nozzle/afterbody geometry, is the topic of this chapter. The relative simplicity of this model 
enabled parametric variations to be implemented and tested easily, some of which are addressed here. 

3.2.1 Experimental apparatus and flow conditions- The experimental model under consideration 
here was analyzed by Cubbage and Kirkham [14] in the late 1960s, and had as its goal the measurement 
of the lift, drag, and pitching-moment increments resulting from the impingement of an underexpanded 
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exhaust plume on a flat plate. The model, pictured in figure 3, consists of an inlet/nozzle assembly attached 
to a flat plate that represents a simplified hypersonic vehicle underbody. The plate is hinged at the nozzle 
exit plane to simulate varying afterbody upsweep angles. The afterbody portion of the model measures 
17.78 cm long by 25.4 cm wide. The nozzle assembly is 9.14 cm wide and 2.29 cm high. Boundary layer 
trips were placed approximately 3 cm downstream of the leading edge to insure turbulent flow over the 
rest of the body. The tests were conducted in the Langley 20-inch Mach 6 tunnel at a stagnation pressure 
of 25 atm and a stagnation temperature of 480 K. The Reynolds number was 20.9x1 0 6 per meter. Pressure 
taps were placed in a grid pattern downstream of the nozzle assembly in the jet expansion region. 

3.2.2 Experimental procedure- Two axisymmetric nozzle geometries were tested with the flat plate 
model, one of which was designed for air (7 = 1 .4), and the other for a gas mixture of 75% sulfur hex- 
aflouride (SF6) and 25% nitrogen (N2) by weight, so chosen to yield an isentropic exponent, 7, of 1.26 at 
the nozzle exit plane. This mixture was chosen because it is relatively inert and stable, producing a dense, 
energetic gas akin to frozen flow in an actual flight condition. The nozzles were designed to produce nearly 
the same momentum flux and Mach number at the nozzle exit plane for both gases, and both had an exit 
diameter of 2.03 cm. A “two-dimensional,” or slot nozzle was also tested with air as the exhaust gas. The 
width of the slot nozzle was the same as the width of the axisymmetric nozzle assembly, 9.14 cm. The 
height of the slot nozzle was 1.52 cm. Air flowing into the inlet of the nozzle assembly was exhausted out- 
side the test section, while the nozzle flow was supplied by a high pressure reservoir. Pressure data were 
reported along the model centerline (between the two inner nozzles) and on the centerline of an outboard 
nozzle. Normal forces on the plate were computed by discrete integration of the pressure data. Measuring 
the power-off static pressure on the plate downstream of the exhaust nozzles determined the power-off 
normal force, with which the researchers were also able to calculate the incremental forces and moments 
on the plate (that is to say, the forces and moments in excess of the power-off values). 

3.3 Numerical Problem Formulation 


Rarely is an actual or proposed problem so simple that it can be modeled exactly in a computational 
analogy. Simplifying the flow physics is essential in light of current computational capabilities; the physical 
approximations made in this study were already discussed in the preceding chapter. In addition, seemingly 
simple geometrical shapes prove astonishingly difficult to model mathematically. In this instance, an oth- 
erwise trivial flat plate geometry was complicated by an assembly of round nozzles, aligned side by side, 
within a square template. To generate a body-conforming grid for such a configuration would have been a 
substantial project in itself. An alternative method that modeled the nozzle assembly with a nonconforming 
grid was found to work quite satisfactorily. This section describes the geometry definition, grid generation, 
and boundary conditions that rendered the numerical analog to the nozzle/afterbody model. 

3.3.1 Geometry definition and grid generation- The parametric variations performed in this exper- 
iment made it possible to assess the effect of afterbody upsweep, nozzle exit pressure ratio, exhaust gas 
isentropic exponent, and nozzle geometry on the forces and moments induced by the impinging exhaust 
plumes. Modeling the same effects computationally presented an opportunity to test the numerical method 
on a very practical application. Geometry definition for the nozzle/afterbody model simply amounted to 
defining a flat plate at arbitrary incidence. The geometry lent itself naturally to a Cartesian grid topology, 
and the three-dimensional grid was generated algebraically. Geometric stretching clustered grid points 
near the solid surfaces. The first grid point away from the plate surface was placed 0.001 cm away, which 
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resulted in an average y + of about 2.5. The quantity y + , a nondimensional distance measured in wall units, 
is defined as 

+ Pw^ru \J Pw t "m y 

y = = ( 3 .2) 

Pw PlV 

where u r is the friction velocity, and r w is the wall shear stress. Wall units make a useful coordinate 
system to analyze turbulent boundary layers, and measuring the grid points in these units provides some 
insight into the grid resolution in terms of the flow variables. In order to make meaningful measurements 
of quantities such as skin friction and heat transfer, it is generally agreed that the normal grid spacing at the 
wall should be less than about ten wall units (private communication from D. S. Chaussee, NASA Ames 
Research Center, Moffett Field,, California). Since the value of y + is not known a priori , experience is 
helpful in choosing appropriate clustering. 

Figure 4 shows a perspective view of the grid topology, which was segmented into two zones. The 
two zones cannot be discerned in the figure because they are one-to-one at the interface. Zone 1 had 32 
points in the stream direction, 33 points across the semispan, and 62 points in the plate-normal direction, 
while zone 2 had 31, 33, and 62, respectively. Excluding the overlapping planes between zones, the grid 
totaled 122,760 points. 

3.3.2 Zonal grid methodology- Dividing the spatial grid into two or more zones offers several advan- 
tages to solving complex flow problems. First, it permits the domain to be segmented into logical subunits, 
such as a forebody, wing, afterbody, and so on. Appropriate placement of zones also eliminates some grid 
singularities, such as discontinuous body cross sections. Finer grids can be used in multiple zone calcula- 
tions than in single zone calculations if the computer’s main memory limits the problem size, because only 
one zone resides in main memory at a time. Finally, the equation set can sometimes be simplified in some 
of the zones, thereby reducing the CPU time required to perform an iteration. 

Several alternatives exist in developing a zonal grid scheme. The simplest is to impose one-to-one 
correspondence of the grid points at the zonal interface, as shown in figure 5. This approach enables direct 
flow field information transfer at each interface grid point, without interpolation whatsoever. While this 
approach is robust and easy to implement, it is not suited to geometrical discontinuities, which occur at 
the nozzle exit plane, for example. On the other hand, the most general zonal scheme permits arbitrary 
overlapping of grids without regard to the interface topology. A scheme with this capability has been 
developed for both two-dimensional and three-dimensional overlaid grids [49]. As yet, the transfer of 
flow field information in this scheme is nonconservative, so shock waves and other flow discontinuities 
may not propagate through zones correctly. Rai has developed and implemented a zonal scheme that 
permits planar zonal interfaces that need not be one-to-one [50]. Furthermore, the scheme ensures that the 
conserved quantities (mass, momentum, and energy) propagate through zones correctly. 

For the present application, zonal grids were introduced primarily to reduce the main memory re- 
quirement. As such, a global grid was first generated about the entire geometry, then it was segmented into 
two zones along the body axis. These zonal interfaces are one-to-one, so no interpolation was necessary 
to pass information from one zone to the next. The zones were overlapped so that three redundant grid 
planes occurred at each interface, enabling boundary conditions at a zonal interface to be obtained from 
the interior of the adjacent zone. 
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3.3.3 Boundary conditions- The solution domain includes only that part of the plate downstream of 
the exhaust nozzles. For the reasons mentioned earlier, the nozzle geometry was not modeled explicitly in 
the grid; rather, the effect of the nozzle flow (and surrounding solid surfaces) was accounted for through the 
inflow boundary conditions, as shown in figure 6. Specifically, a logical array was defined for the inflow 
plane that determined the type of boundary condition to apply at each point in that plane. Flow variables 
at points outside of the nozzle assembly region were fixed from a previously calculated turbulent flat plate 
solution. Flow variables at points within the nozzle assembly region of the inflow plane were specified to 
be either: 1) nozzle exhaust conditions, with all conservative variables fixed, or 2) solid surface conditions, 
with no slip and dT/dn = 0 , dp/dn = 0 . The temperature boundary condition, dT/dn = 0 , corresponds 
to an adiabatic wall condition, and is implemented in the code by specifying T\ - Tz- The pressure 
boundary condition corresponds to the normal momentum equation at the wall, and is implemented by 
specifying pi = P 2 . Both of these conditions are strictly correct only for a grid that is orthogonal at 
the wall. Having fixed the temperature and pressure at the wall, the conservative variables p and e were 
decoded using the equation of state. 

Since the geometry and flow field are symmetric about the centerline of the model, only one half of the 
model was included in the domain. Symmetry boundary conditions were applied at the plate centerline: all 
conservative variables were reflected symmetrically acrosss the symmetry plane, except the y-momentum 
(pv), which was reflected antisymetrically. Because the Osher scheme uses a five-point difference stencil, 
two reflection planes are required at the symmetry plane. Far field conditions were specified on the outer 
boundaries, in which all flow variables are set to their freestream values. At the outflow boundary, the con- 
servative variables were linearly extrapolated from the points upstream. No-slip conditions were imposed 
at the plate surface, and turbulent flow was assumed from the inflow plane. 

3.4 Solution Procedure 


The conservative variables at interior grid points must be initialized to start the iterative solution 
algorithm. When no better initial guess was available, the flow variables were initialized to their freestream 
quantities. Since this “guess” of the solution is far from the converged result, a large perturbation occurs at 
the first time step. Stability bounds constrained the time step to a small value initially. It was then gradually 
increased as the solution progressed to the steady-state answer. Rather than fixing the time step explicitly, 
the UWIN code allows specification of the Courant number, defined as 

CNBR = r— — (3.3) 

^ max 

where \ max is the largest eigenvalue of the flux Jacobian matrices. The time step is then calculated from the 
prescribed Courant number and the maximum eigenvalue. This method is convenient because the Courant 
number inherently accounts for the grid and flow features that affect stability. Based on experience with 
the UWIN code, the Courant number was initially set to a value of 0.1. As the solution progressed, the 
Courant number was gradually increased to a maximum value of 100. Although the code demonstrated 
stability at Courant numbers as high as 500, generally no benefit to the convergence rate was observed at 
values greater than 100. 

Two criteria determined convergence: 1) the Lj norm of the residual was reduced three orders of 
magnitude from its maximum value, and 2) the maximum change in density in one time step was less 
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than 10 -4 . Upon satisfying both of these criteria, the Courant number was reduced to a value of unity 
and the solution was iterated further to eliminate high frequency errors. In practice, the axisymmetric 
nozzle cases did not satisfy the second convergence criterion. This is thought to be due to unsteadiness 
caused by periodic vortex shedding from the nozzle lip. Evaluation of the quasi-steady results indicate no 
demonstrable effect of the unsteadiness on the results presented here. 

Solutions started from freestream conditions generally converged in about 600 iterations. Solutions 
were run on both the Cray X-MP and Cray 2; the code expends about 144 microseconds of CPU time 
per iteration per grid point on either. This equates to approximately 50 megaflops (million floating point 
operations per second) for the UWIN code. Total CPU time for a solution was about four hours. 


To expedite convergence, in many instances new solutions were generated by perturbing converged 
solutions. For example, to calculate the flow solution at a new exit pressure ratio, the inflow boundary con- 
ditions were modified to reflect the new jet exhaust variables, and flow variables at the interior points were 
initialized from the previous solution. This procedure reduced run time by about 50%, thereby enabling a 
more comprehensive analysis of parametric variations. 

3.5 The Nozzle/Afterbody Flow Field 


To facilitate later discussions, it is helpful first to describe the important features of an underexpanded 
jet flow, along with the peculiarities of the nozzle/afterbody flow field. A jet is termed underexpanded 
when the static pressure at the nozzle exit plane is greater than the ambient static pressure. Upon leaving 
the nozzle, the flow “expands” by turning toward the region of lower pressure, as shown in figure 7. The 
adjacent supersonic external flow is forced to turn also, and an oblique shock wave (the plume shock) 
emanates from the nozzle lip. Meanwhile, an expansion fan propagates into the plume flow, eventually 
contacting and reflecting off of the afterbody surface. For two-dimensional and axisymmetric flows, waves 
continue to reflect off of the contact surface and nozzle centerline, giving rise to the diamond pattern 
sometimes seen in rocket plumes. For nonaxisymmetric geometries, waves propagate out of phase with one 
another. The ensuing nonsimple wave interactions dissipate the pressure gradients through fluid viscosity, 
so the wave structure is less apparent. The flow is yet more complex for the nozzle/afterbody model, as 
shown schematically in figure 8, owing to the side-by-side axisymmetric nozzle configuration. The flow 
stagnates between adjacent nozzles. Expansion waves propagate into the circular jet flow from the entire 
circumference of the nozzle. And when adjacent nozzle plumes interfere with each other, new shocks form. 

This flow pattern is also evident in the numerical solutions. For example, figure 9 shows the calculated 
pressure contours on the plate, at the symmetry plane of the model, and at the inflow plane. The freestream 
Mach number, Moo * for this solution is six, as for all the other nozzle/afterbody solutions presented. The 
angle of attack, a, is five degrees, and the Reynolds number. Re, is 209,000 cm -1 . The afterbody upsweep 
angle, designated by e, is zero degrees. The nozzle exit pressure ratio (p e /Poo) is 6.23, and the exit Mach 
number M e is 3.5. At the inflow plane, pressure contours outline the nozzle circumferences. As described 
earlier, a Cartesian grid topology is employed with explicit inflow boundary conditions. The discrete 
representation of the round nozzles is evident from the pattern of pressure contours at the inflow plane. 
The plume shock is also readily evident in this view at the symmetry plane. Pressure peaks on the plate 
are apparent at the nozzle lip (where the exhaust flow first impinges on the plate), and further downstream 
where adjacent plumes first interact. 
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Closer inspection of the flow near the nozzle exit reveals large regions of recirculating flow. A detailed 
view of the limiting streamlines near the plate for this case is shown in figure 10. By restricting these 
particles to remain in the coordinate surface nearest the plate, the path they trace is akin to an experimental 
oil flow visualization [51]. Between the nozzles, the flow stagnates for a distance of approximately one 
nozzle diameter downstream of the exit plane. At the downstream edge of this region, there exists a saddle 
point singularity where fluid particles are either entrained by the plume flow and carried downstream, or 
drawn into the eddying flow between the nozzles. The flow in the stagnant region then separates along a pair 
of separation lines, both of which terminate at the solid surface of the nozzle template. The postulated skin 
friction pattern for this flow is included in figure 10. According to rules of skin friction topology, the number 
of nodal points must exceed the number of saddle points by exactly two. As the figure indicates, the model 
features nine nodal points of attachment (one for each nozzle, one on each side of the nozzle assembly, 
and one in each aspiration zone), and one nodal point of separation (at the downstream boundary). Also 
shown are eight saddle points that result from the recirculating flow near the nozzle exit plane. Thus, the 
topological rule is satisfied and the interpretation is physically possible. However, closer analysis of this 
region may reveal yet more flow details. The pressure in the stagnant region is close to the ambient static 
pressure, and hence no lift or thrust is obtained here. 

Surface pressure coefficients are compared with the experimental data for this case in figure 11. In 
figure 11(a), the pressure along the model centerline is seen first to rise due to a shock at x/d pa 1 .5, 
then fall to near ambient pressure at x/d pa 5. The shock is the result of adjacent plumes interfering 
with one another as they expand to ambient pressure. The flow then overexpands, and at x/d pa 6, a 
recompression is apparent The comparison in this figure is generally good, though the strength of the 
first shock is underpredicted somewhat. This is due in part to the discrete modeling of the round nozzles, 
as will be shown later in a study of grid refinement: more accurate modeling of the nozzle results in 
greater interference and a higher pressure peak. Figure 11(b) shows the surface pressure coefficients on 
the centerline of the outboard nozzle. The pressure falls rapidly from its value at the exit plane, then 
reaches a plateau which occurs at x/d pa 1 in the experimental data, and at x/d pa 2 in the calculated 
result. The pressure plateau occurs where the expansion of the exhaust plume is restricted by interference 
with the adjacent plume. Downstream of this point, the flow becomes slightly overexpanded, then recovers 
to ambient pressure. Discrepancies are again attributed to grid resolution at the nozzle exit plane. 

In the discussion above, reference was made to the high pressure region on the plate downstream of 
the nozzles. The term “footprint” is commonly used to describe the extent to which higher than ambient 
pressure prevails. As can be seen in figure 9, most of the normal force on the plate comes from the pressure 
peaks associated with the initial impingement of the exhaust fluid on the plate and, further downstream, 
the pressure peaks associated with mutual plume interference. These pressure peaks exhibit part of the net 
force-producing footprint on the afterbody. A useful parameter with which to compare afterbody forces is 
the normal force coefficient, designated c n and defined as 




where N is the force normal to the afterbody, and A re f is the reference area of the plate (452 cm 2 ). For 
this case, the normal force coefficient on the plate was computed to be 0.1256. Since a = 5° and € = 0° 
for this case, the afterbody force coefficients in Cartesian components are c x = 0 .0109 and c z = 0 . 125 1 . 
The quantity c n will be useful later in correlating the calculations with experimental results. 

3.6 The Influence of Parametric Variations on Afterbody Forces 


Having gained some appreciation for the dynamics of underexpanded jet impingement, it is now of 
interest to investigate parametric variations of nozzle operating conditions. As parameters such as the exit 
pressure ratio, afterbody upsweep, and nozzle flow conditions change, the plume/afterbody interaction will 
be altered, impacting overall performance. Becoming familiar with these performance derivatives will be 
useful in interpreting later results. 

3.6.1 Effect of nozzle exit pressure- Whereas turbojet engines derive most of their thrust by ex- 
panding the captured air to very high velocity inside the nozzle, scramjet engines must rely on external 
expansion of the fluid. A scramjet nozzle that completely expanded the exhaust gas could feasibly be 
larger than the vehicle it was designed to propel! Thus, fluid leaving the scramjet nozzle will be highly 
underexpanded. The nozzle exit pressure ratio — defined here as the ratio of static pressure at the nozzle 
exit plane to the freestream static pressure — is an important determinant of the afterbody forces. It influ- 
ences the shape and force of the footprint the exhaust plume makes as it impinges on the afterbody, thereby 
contributing to the overall lift, thrust, and moment 


Solutions for the nozzle/afterbody model were obtained at an otherwise fixed operating condition 
with variations in exit pressure ratio. Comparisons of the surface pressure coefficient at two values of exit 
pressure ratio are shown in figure 12 for a = 5°, e = 0°, and p e /p<x> = 2 .07 and 6 .23. The calculated 
results show good agreement with the data except in predicting the location of plume interference as noted 
earlier. Note that the p e /p«> = 1 data — the condition at which the nozzle flow is completely expanded 
— indicates that the flow initially overexpands because of low pressure in the aspiration region near the 
nozzle exit plane. Further downstream, recompression to ambient pressure is observed. 


To assess the sensitivity of the afterbody normal force to the exit pressure ratio, solutions were ob- 
tained for p e /poo = 1,2,3, and 4 at a = 0°, c = 0°. Figure 13 presents incremental normal force 
coefficients as a function of the nozzle exit pressure ratio. To compute the “incremental” value, the experi- 
mentalists subtracted the power-off normal force coefficient, c,*, , from the power-on value. The numerical 
solutions were correlated with this data by normalizing the calculations about p e /poo = 1 , thereby deducing 
c,*, as follows: 


M 


numerical 
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(3.5) 


The linear relationship between exit pressure and normal force is tentatively explained as follows: 
an initial expansion wave propagates into the plume with a slope determined only by the plume Mach 
number. In the region of exhaust flow undisturbed by expansion waves, the pressure is nearly the same as 
the exit pressure; thus, pressure on the plate will persist at this value until the first expansion wave reaches 
the surface. Thereafter, expansion proceeds until the ambient pressure is reached. Since the exit Mach 
number was held constant for these solutions, the expansion began at approximately the same location for 
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all cases. Thus, as the exit pressure ratio was increased, the pressure in the basic footprint region rose 
proportionately, but the area over which it acted remained largely unchanged The result was a linear 
increase in the normal force coefficient with exit pressure ratio. 

This hypothesis is verified by comparisons of the footprint shape at each pressure ratio, shown in figure 
14. The footprint shape is defined for this illustration by the location of the pressure contour on the plate 
that is 10% above freestream pressure. Hence, this contour envelops the region on the plate that produces 
force in excess of power-off conditions. At p e /poo = 3 and 4, the contours are very comparable, showing 
that the exit pressure has only a small influence on the size and shape of the footprint. At p e /p<x> = 2 , the 
footprint is considerably smaller, and at p e /poo = 1 , the footprint is virtually nil. However, these two lower 
values are so close to the ambient pressure even at the exit plane, that small pressure variations (arising 
from plume interference, for example) affect the footprint shape substantially. 

3.6.2 Effect of afterbody upsweep- The increase in thrust derived from an upswept afterbody has 
already been discussed, though no mention was made of how much upsweep optimizes a design. The 
schematic in figure 15 shows that, on one hand, upsweep increases the overall footprint area of the exhaust 
plume. However, the flow expands in turning to follow the afterbody contour, so the footprint area on an 
upswept afterbody will exhibit lower pressure than the unswept case. Too much upsweep results in an 
overexpanded flow that is likely to separate from the body, reducing lift and increasing drag. Clearly, the 
optimal upsweep angle depends closely on the design exit pressure ratio. 

Consider the surface pressure contour comparisons at a = 0°, p e /poo = 3 shown in figure 16. The 
footprint on the plate at e = 4° extends to x/d « 5, whereas at e = 0°, it reaches only to x/d & 4. A two- 
dimensional method of characteristics analysis predicted values of 5.1 and 4.2, respectively. Furthermore, 
because the upswept afterbody permits the exhaust flow to expand partially, the interference peaks for 
c = 4 0 are not as high as the e = 0 0 case. Table 1 below compares the forces exerted on the plate for 
this case and for p e /Poo = 1 • One interesting observation from the table is that the incremental normal 
force is nearly the same for the flat and upswept afterbodies at a given exit pressure ratio. Apparently, the 
phenomena associated with upsweep are self-cancelling in these cases. The upswept afterbody’s larger 
footprint tends to add to the plate force, while the added expansion decreases the pressure of the exhaust 
gas and tends to diminish the plate force. However, since the external flow must turn also, the ambient 
pressure drops along with the exhaust pressure. In fact, since the external flow is at Moo = 6 , Prandtl- 
Meyer theory predicts that the external flow turning four degrees will experience a 59% drop in pressure; 
the exhaust flow, which is at M e = 3 .5 , will drop by only 36% [52]. The results in Table 1 suggest that 
the net effect is a slight loss in incremental normal force. 


Table 1. Force coefficients for the flat and upswept afterbody. Moo = 6 , a = 0°. 


Pe/ Poo 

l 

3 

€ 

0° 

4° 

0° 

4° 

C* 

.0000 

.0013 

.0000 

.0021 

Cz 

.0225 

.0190 

.0338 

.0296 

c n 

.0225 

.0190 

.0338 

.0297 
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3.6.3 Effect of exhaust gas composition- At high altitude and Mach number flight conditions, scram- 
jet exhaust species will continue to undergo chemical reactions after leaving the combustion chamber. 
However, if expansion occurs too rapidly, the temperature may drop below the threshold necessary for 
recombination reactions to proceed, freezing the composition in an intermediate state (private communi- 
cation from G. S. Deiwert, NASA Ames Research Center, Moffett Field, California). When this happens, 
the exhaust flow traverses the body without reacting, behaving as an ideal gas that mixes with another 
ideal gas, the freestream air. By coupling a species continuity equation to the Navier-Stokes equations as 
described in chapter 2 , it was possible to assess the effect that different (ideal) exhaust gases have on the 
afterbody forces. While there is no simplified theory to analyze such a flow, a heuristic argument is put 
forth here. 

The experiment to which the present results are compared tested a mixture of 75% sulfur hexaflouride 
(SF 6 ) and 25% nitrogen (N 2 ), for which 7 e = 1 .26 , exhausting into air (7 = 1 .4 ). Kinetic theory [46] 
explains sulfur hexaflouride ’s lower value of 7 as the result of its more numerous modes of internal energy 
than those of air. At moderate temperatures, air, composed almost entirely of diatomic molecules, exhibits 
five modes of internal energy: three translational modes, and two rotational modes. The kinetic model 
predicts that 7 = ( n+ 2) /n, where n is the number of internal energy modes. Thus, 7 air = 1 .4 , as verified 
by observation. However, SF6 is a much larger molecule possessing numerous vibrational modes, many 
of which are fully excited at room temperature. This results in a lower value of 7. Internal energy must 
be distributed over many more modes in SFe than air; therefore, SF6 has more internal energy than air 
at the same temperature. As the gas expands and cools, energy cascades out of the other internal energy 
modes and into kinetic energy. This will reveal itself in a smaller drop in pressure for a given volumetric 
expansion. Conversely, expanding from a fixed exit pressure to some lower pressure will require a greater 
volumetric expansion of 7 e = 1 .26 fluid than of 7 e = 1 .4 fluid. Hence, the plume footprint is expected to 
be larger, and the expansion more gradual for SF6 than for air. 

Figure 17 bears out this argument. At a = 0 °, c = 0°, p e /poo = 4 , the footprint for SF 6 exhaust gas 
is notably larger than for air: it persists an additional nozzle diameter downstream, and it spreads laterally 
by about another one-half nozzle diameter. Figure 18 shows pressure contours on the plate, demonstrating 
that the pressure within the SF 6 footprint is higher as well; for this case, the incremental normal force 
for 7 e = 1 .26 was about 25% higher than for 7 e = 1 .4. Comparisons of the calculated and measured 
normal force coefficient at several exit pressure ratios are shown in figure 19. The variation of normal 
force coefficient with exit pressure ratio is still linear, but the slope is about 30% greater for 7 *. = 1 .26 
than for 7 e = 1.4. Therefore, at the same exit conditions, a gas with a low value of 7 will produce greater 
afterbody forces than one with a high value of 7 . 

Contours of the exhaust gas species concentration are shown in figure 20. Since convection is the 
only transport mechanism modeled, the contours retain the general shape of the contours at the nozzle 
exit plane. They distort somewhat due to the expansion of the plumes, but the most significant species 
transport is seen to occur near the plate. In the boundary layer near the surface, lateral spreading is markedly 
enhanced because the velocity is so much smaller than the flow outside the boundary layer. The crossflow 
velocity becomes the same order of magnitude as the stream wise velocity, and exhaust gas species are thus 
convected much further laterally in the boundary layer than outside of it. 

3.6.4 Effect of nozzle geometry- In figure 10, it was shown that adjacent axisymmetric nozzles leave 
regions of stagnant flow near the exit plane. This observation begs the question, “does a nozzle geometry 
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that produces a ‘cleaner’ flow field near the exit generate greater afterbody forces?” To investigate this 
question, an alternate “slot” nozzle was tested against the axisymmetric nozzle geometry. The dimensions 
of the slot nozzle model were determined by specifying the overall width of the nozzle assembly to remain 
the same. The nozzle height was then fixed so as to produce the same momentum flux through the slot 
nozzle as the four axisymmetric nozzles at the same exit conditions. Test conditions for the comparison 
were a = 5 e = 0 °, and p e /poo = 6 .23 . 

Comparisons of the calculated surface pressure coefficients with experimental data are shown in figure 
21. Centerline pressure for the slot nozzle is seen to be consistently underpredicted from x/d = 0 to 
x/d & 5. This is thought to be a result of two discrepancies between the experiment and the numerical 
model: 1) nonuniform experimental exit flow conditions, causing a pressure peak near the exit face, and 
2) geometrical disparity, causing the predicted onset of expansion to occur about one nozzle diameter 
upstream of the experimentally determined location. 

The flow pattern near the plate renders a clearer visualization of why the two nozzle geometries per- 
form differently. Figure 22 presents limiting streamlines for the two nozzle geometries. These streamlines 
are colored according to the local static pressure, adding insight as to how the two flow fields generate 
plate forces. The slot nozzle produces a virtually two-dimensional flowfield away from the edges of the 
nozzle assembly. The flow is, indeed, much “cleaner” than the axisymmetric geometry, showing no re- 
circulating flow near the nozzle exit. The payoff is in the afterbody normal force: Cn^ = 0 .0778 , while 
Ci^xx ~ 0 0409 . The recirculation regions on the axisymmetric geometry dissipate nearly 50% of the in- 
cremental normal force compared to the slot nozzle! As suggested by figure 22, most of the additional 
force is generated where the axisymmetric nozzle flow creates low pressure recirculation regions. 

3.7 The Effect of Grid Resolution on Computed Results 


As a necessary final step to the nozzle/afterbody investigation, the convergence of the numerical 
results is tested on a series of grids. Any finite-difference method carries with it truncation errors arising 
from approximations to the differentiations. Furthermore, it can be shown that the truncation error is 
proportional to the cell volume (lecture notes from H. Lomax, NASA Ames Research Center, Moffett 
Field, California). Hence, as the grid is refined, the cell volume decreases, and truncation errors diminish. 
Obtaining numerical solutions on a series of grids tests the accuracy and uniqueness of the results. The 
value of some relevant flow field variable should approach a limit asymptotically as the grid is refined. 
This section presents the results of a grid refinement study performed for the nozzle/afterbody model at 
a = 5°, e = 0°, and pe/poo = 6 .23. 

To minimize the effects of unequal grid stretching on the results of this study, results were first obtained 
on a fine grid, which was then coarsened twice by discarding alternate grid points in each of the three 
coordinate directions. Though this technique does not preserve identical stretching on the coarsened grid, 
it does produce proportionately similar grid cells from one grid to the next. Discarding alternate grid cells 
in all three directions results in an eight-fold reduction in the total number of grid points. The fine grid 
dimensions were 119 points in the streamwise direction, 63 in the spanwise direction, and 123 in the plate- 
normal direction, for an overall total of 922 ,131 points. Hence, the medium grid dimensions were 60, 33, 
and 62 points, respectively, for a total of 122 , 760 points; the coarse grid had 31, 17, and 31 for 16 , 337 . 
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Pressure coefficient comparisons among solutions on the three grids are shown in figure 23. The 
coarse grid result follows the broad trends in the data, but fails to place the interference peak correctly in 
figure 23(a), and shows no evidence of recompression at x/d = 6 . The medium and fine grid results are 
nearly identical to one another, the only noteworthy difference being that the fine grid result shows a higher 
interference peak (figure 23(a)), with a curious saddle profile. Figure 23(b) shows similar trends. 

Figure 24 compares boundary layer profiles among the solutions. The profiles between adjacent noz- 
zles, figure 24(a), shows some differences in the solutions. The coarse grid predicts attached flow at this 
station, while the medium and fine grids show significant reverse flow. The flow reversal is precipitated 
by shock waves from the interaction of adjacent plumes, and the discrepancies seen in this figure are most 
likely a result of how well the plume interaction process is resolved. On the nozzle centerline, figure 24(b), 
all the results are identical, indicating that even the coarse grid adequately models the attached flow. To 
resolve the viscosity-dominated flow field near the nozzle exit would likely require a grid much finer than 
those studied here. 

However, plate force measurements would not be sensitive to further refinement of the grid. Plate 
normal force coefficients, shown in Table 2, indicate an asymptotic trend toward a value of about 0.138. 
The medium grid dimensions, which were used for the other nozzle/afterbody grids as well, appear to 
resolve most of the flow characteristics adequately, and can be used to calculate integrated quantities to 
within about 10% of the asymptotic value. 


Table 2. Effect of grid density on plate normal force. 
Moo = 6 , at =5°, e = 0°, p e /poo = 6 .23. 


Grid 

Dimensions 

Total Points 

Cn 

Coarse 

h ““ 31 - 17 - 31 

16,337 

0.1009 

Medium 

60 • 33 • 62 

122,760 

0.1256 

Fine 

119 • 63 • 123 

922,131 

0.1362 


3.8 Summary 


The nozzle/afterbody model has provided a reasonable idealization of the flow about an actual vehicle 
afterbody. By restricting the domain to only the region of interest, computational requirements were re- 
duced enough to permit solutions for several parametric variations. Numerical results compared favorably 
with the experimental data, and inadequate resolution of the nozzle geometry may be the principal source 
of discrepancies. Having calibrated the code and gained an understanding of the exhaust plume/afterbody 
flow field, attention is now turned to the effect of the exhaust plume on a realistic hypersonic vehicle 
configuration. 
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Figure 1. Schematic representation of the scramjet propulsion cycle. 
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Figure 2. Control volume of nozzle and afterbody for force calculations. Si is 
the nozzle exit face, S* is the afterbody surface. 
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(a) front perspective 


NOZZLES 



(b) rear perspective 

Figure 3. Nozzle/afterbody model geometry. 
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Figure 4. Cartesian grid topology for nozzle/afterbody model. 
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INFLOW BOUNDARY CONDITION 
O TURBULENT FLAT PLATE SOLUTION 
A SOLID SURFACE BOUNDARY CONDITIONS 
O JET EXHAUST CONDITIONS 
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Figure 6. Inflow boundary conditions for nozzle/afterbody model. 
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Figure 7. Schematic of an underexpanded jet flow field. 
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NOZZLE EXTERNAL FLOW 



(a) overall view 



(b) detail near exit 

Figure 8. Schematic of the interacting axisymmetric nozzle flow near plate. 
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Figure 9. Perspective view of pressure contours on plate surface, at the symme- 
try plane, and at the inflow plane. =6, a = 5°, 

e = 0°, Pel Poo = 3. 
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(b) postulated skin fnction lines 

Figure 10. Detail of surface flow pattern near the nozzle exit plane. = 6, 
a = 5°, e = 0°, Pr/poo = 3. 
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(a) between nozzles 



(b) nozzle centerline 


Figure 11. Comparison 
cients. Moo 


of calculated and experimental surface pressure coeffi 
= 6, a = 5 , e = 0 , Pe/Poo ~ 3. 
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(ft) between nozzles 
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x/d 


(b) nozzle centerline 

Figure 12. Variation of surface pressure coefficient with exit pressure ratio. 
Moo = 6, a = 0°, e = 0° 
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Figure 14. Comparison of footprint shapes for a progression of exit pressure ra- 
tios. Moo = 6, a = 0°, e = 0°. 
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FOOTPRINT REGION 



(a) without upsweep 



(b) with upsweep 

Figure 15. Schematic comparison of an underexpanded nozzle flow field with and 
without upswept afterbody. 
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Figure 16. Comparison of surface pressure contours for flat and upswept after- 
body profiles. Moo =6, a = 0°, p e /poo = 3. 
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Figure 17. Comparison of footprint shapes for different exhaust gases. M «, = 6, 
a = 0°,c = 0°, Pe/Poo = 4. 
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Figure 18. Comparison of surface pressure contours for different exhaust gases. 
Moo = 6, a = 0°, e = 0°, p e /poo = 4. 



Figure 19. Perspective view of exhaust gas concentration contours. M m = 6, 
a = 0°, £ = 0°, p e /poo = 4, 7, = 1.26. 
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(a) axisymmetric nozzles 
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x/d 


(b) slot nozzle 

Figure 21. Comparison of calculated and experimental surface pressure coeffi- 
cients for axisymmetric and slot nozzle geometries. Moo — 6, 
a = 5°, e = 0°. 
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Figure 22. Comparison of limiting streamlines colored by pressure near noz- 
zle exit for axisymmetric and slot nozzle geometries. Afoo = 6, 
a — 5°, e = 0°, Pt/Poo = 6.23. 
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MEDIUM GRID COMPUTATION 
FINE GRID COMPUTATION 



(a) between nozzles 



(b) nozzle centerline 

Figure 23. Comparison of calculated and experimental surface pressure coeffi 
dents for grid refinement study. =6, a = 5°, e = 0°, 

Pe/Poo = 6.23. 
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(a) between nozzles 
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(b) nozzle centerline 

Figure 24. Comparison of boundary layer profiles at x/d = 1 for grid refinement 
study. Moo =6, a = 5°, t = 0°, Pe/Poo — 6.23. 
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CHAPTER FOUR 


ANALYSIS OF GENERIC HYPERSONIC VEHICLE FLOW FIELDS 


4.1 Introduction 


In the preceding chapter, it was shown that propulsion system design parameters have a large influence 
on the afterbody flow field and vehicle forces. It is now appropriate to address questions relating to installed 
scramjet performance: how does the impinging exhaust plume affect vehicle balance? Can proposed nozzle 
exit conditions produce enough thrust to overcome drag? How does the exhaust plume interact with a 
nonuniform external flow? What impact does thrust vectoring have on vehicle forces and moments? These 
are some of the questions that will be addressed in this chapter. Investigating these issues required a realistic 
hypersonic vehicle geometry, including an exhaust nozzle, with experimental data available for a power-on 
flight condition. Though experimental activity in this area is in progress (see reference [20], for example), 
little data already exists. The fact that most candidate experiments are classified exacerbated the difficulty 
in locating data. The generic hypersonic vehicle from Lockheed-Georgia (private communication with D. 
Schuster, Lockheed-Georgia Co., Marietta, Georgia)) best filled these requirements. 

The generic hypersonic vehicle geometry provided a unique opportunity to solve power-on flight 
conditions for a realistic configuration, and to compare those results with experimental data. This chapter 
presents the results of power-on solutions for this body at supersonic and hypersonic speeds. The exper- 
iment included tests at transonic and supersonic conditions; experimental apparatus and procedures for 
these tests are described first. A discussion of formulating the numerical problem includes surface geom- 
etry definition, grid generation, and boundary conditions. Three-dimensional grid generation by means of 
a two-dimensional grid generator is summarized. Also, the boundary conditions and solution procedure 
presented in the previous chapter are modified to accomodate this geometry as explained below. 

Flow field analyses for the generic hypersonic vehicle are presented in two segments. Supersonic 
(Moo = 1 .46 ) solutions at a series of exit pressure ratios are compared with experimental data. The 
afterbody flow field is described, and the effect of the exhaust plume on vehicle forces and moments is 
enumerated. Next, solutions are presented at a practical hypersonic cruise condition (Moo = 6). These 
results demonstrate that the flow field of this strongly interacting exhaust plume is rich in information about 
afterbody aerodynamics. The chapter concludes with a study of thrust vectoring as a possible means to 
balance the vehicle. 


4.2 Experimental Apparatus and Procedure 


The generic hypersonic vehicle model was designed and tested by Lockheed-Georgia to study the 
transonic and supersonic performance of a NASP-like configuration. The model, pictured in figure 25, 
features a pointed nose, blended wing-body, and exhaust nozzle with afterbody contouring. The overall 
vehicle length is 50 cm, and the wing semi-span reaches a maximum of 10 cm. In order to calculate force 
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and moment coeficients for the numerical solutions, a reference area and reference length were required. 
The approximate planform area (500 cm 2 ) and body length (50 cm) provided these values. The afterbody 
upsweep at the centerline is 20 degrees, and the nozzle is 0.5 cm high. Power-on simulations were ac- 
complished by blowing air through the nozzle. The model was supported by a blade mount attached to 
the upper surface; the wing shielded disturbances created by the mount so that it wquld not interfere with 
the nozzle flow. High pressure lines for the nozzle effluent were ducted through the mount to an internal 
plenum. Pressure taps were placed on the model lower surface downstream of the nozzle, on the centerline 
and at a few off-centerline locations. 

The tests were performed in the Lockheed 20 by 28-inch compressible flow tunnel at ffeestream Mach 
numbers from 0.8 to 1.46. At = 1.46, tests were conducted for nozzle exit pressure ratios from 0 
(power off) to 6.27. The tests employed a “cold” jet, meaning that its stagnation temperature was equal to 
the freestream stagnation temperature. By design, the exhaust flow was to be choked at the exit plane, but 
manufacturing imperfections caused the flow to choke inside the model. Exit plane measurements showed 
that the exit Mach number varied circumferentially, and ranged approximately from 1 . 1 2 to 1 .2. The tunnel 
stagnation pressure was 1.4 atm, stagnation temperature was 396 K, and the Reynolds number was 230,000 
cm -1 . 


4.3 Numerical Problem Formulation 


4.3.1 Geometry definition and grid generation- The geometry definition for this model was sup- 
plied in the form of body-axis-normal cross section coordinates at specified stations along the length of 
the body. To produce a smooth grid with adequate resolution, it was necessary to generate new cross sec- 
tions by interpolating between succesive body stations. Not only did this yield a smoother distribution of 
grid points in the streamwise direction, but it also improved resolution of the afterbody downstream of the 
nozzle, where rapid plume expansion occurs. 

The three-dimensional grid was generated using a two-dimensional hyperbolic grid generation pro- 
cedure [53] at each cross section. This technique solves a system of hyperbolic equations for the grid 
coordinates (y, z) in the cross sectional plane. Because the system of equations is hyperbolic, the solution 
procedure marches from the body outward to the outer boundary. Consequently, generating grids is com- 
putationally efficient (approximately 50 /xsec of Cray XMP CPU time per grid point for these grids), and 
improvements or alterations to the grid are accomplished quickly. 

The two-dimensional grids were then assembled to form a three-dimensional grid about the entire 
body, producing a cylindrical, or O-H, grid topology as shown in figure 26. This method and grid topology 
are quite amenable to blended bodies. At the nozzle exit plane, where the geometry has two cross sections, 
grids were generated for the section both including and excluding the nozzle contour, as shown in figure 
27. This resulted in a planar zonal interface that was not one-to-one. The UWIN code features a trilinear 
interpolation routine wherein coefficients are computed at each interface grid point. Since the grids are 
stationary, these coefficients are computed only once, and are subsequently used to transfer information 
across the interface. 

The grid was then segmented into four sequential zones, each of which had 16 points in the stream 
direction, 89 points in the circumferential direction, and 50 points in the radial direction. Disregarding 
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the overlap planes, the computational domain consisted of 244,750 grid points. Exponential clustering 
in the radial direction resulted in a normal spacing of the first point off the body of about 0.01 cm. At 
Moo = 1 .46 , this spacing equated to y + « 5 , while at Moo = 6, j/ + was about 1. 

4.3.2 Boundary Conditions- For bodies with sharp or blunt noses in supersonic flow, boundary 
condition procedures customarily require that the upstream boundary be placed sufficiently far from the 
body to capture the bow shock. However, the formulation of the Osher scheme in the UWIN code possesses 
a known problem (private communication with J. Flores, NASA Ames Research Center, Moffett Field, 
California) in solving the flow along the singular grid line upstream of the body. Work on this problem is 
in progress, but in the meantime an alternative is needed. One possibility is to generate a starting solution 
with another code, and use that as fixed inflow data to the UWIN code at a body station downstream of 
the nose. This method is cumbersome, though, and since codes make different approximations and errors, 
information transfer between codes introduces errors at the inflow plane. 

A viable alternative for this study was devised by placing the upstream boundary at the pointed nose 
of the body. Thus, the inflow plane was a true polar grid, with the only point on the body being at the 
origin. Flow conditions on this plane were fixed to be ffeestream at every point except the origin, where 
flow conditions were specified based on tabulated normal shock conditions [52]. These conditions make 
sense from a physical standpoint because in supersonic flow, disturbances propagate only downstream; 
hence at the nose, the only point where the flow “knows” about the body is the origin. Evaluation of the 
results indicates that this boundary condition did indeed produce a realistic inflow condition. 

The other boundary conditions were specified as before: symmetry conditions at the windward and 
leeward pitch planes of symmetry, ffeestream at the outer boundary, no slip and dp/dn = 0 , dT/dn = 0 
at the body, and supersonic outflow at the downstream boundary. At the nozzle exit plane, flow variables 
at points within the nozzle circumference were specified from the exhaust conditions, while variables at 
points outside the nozzle were specified by interpolation from the interface plane of the upstream zone. 

4.4 Solution Procedure 


While supersonic flows over streamlined bodies such as this can be solved efficiently by parabolized 
Navier-Stokes methods [31], these codes cannot solve flows with streamwise separation. As will be shown 
shortly, interactions between the exhaust plume and the external flow do precipitate streamwise separation 
in some instances. Also, part of the intent of this work was to demonstrate a single code, tip-to-tail solution 
capability. As such, the UWIN code solved the entire flow field around the generic hypersonic vehicle. 
First, a power-off solution was generated by imposing no-slip conditions on the face of the exhaust nozzle. 
Transition was assumed at the nose, with turbulent flow prevailing over the entire body surface. Subse- 
quent power-on solutions were obtained by specifying nozzle exhaust conditions, and iterating the two 
downstream zones. Iterating the two upstream zones was unnecessary because in these flow fields, up- 
stream influence of the nozzle is localized. The effect of the exhaust plume on the upstream flow was 
easily contained within the zone immediately upstream of the nozzle. 

Solutions were run on the Cray X-MP, and the initial power-off solutions converged in about 400 iter- 
ations, requiring about 4 hours of CPU time. Subsequent power-on solutions required about 500 iterations 
of just the two downstream zones, taking about 2.8 CPU hours. 
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4.5 Generic Hypersonic Vehicle Flow Solutions at Moo = 1 .46 


In order to correlate the numerical results with available experimental data, the generic hypersonic 
vehicle geometry was first studied at a freestream Mach number of 1.46 and an angle of attack, a, of zero 
degrees. A series of exit pressure ratios were tested, from p e /Pc© = 0 (power off) to p e /poo = 6 .27. 
This section presents a description of the flow fields, comparisons to experimental data, and determination 
of vehicle forces. To begin, pressure contours for the entire solution domain are shown in figure 28 at 
Pe/Poo = 6 .27. Many interesting phenomena apparent in this flow field merit further study: formation of 
shocks at the bow and at the wing, and their eventual interaction; heat transfer peaks at the nose and sharp 
leading edge of the wing; and the flow pattern around the square wing tip, to name a few. Unfortunately, 
these aspects are beyond the scope of this investigation, and will be set aside. Instead, the windward-side 
flow field on the afterbody is of current interest, and the details of this region are now explored. 

Figure 29 shows a detailed view of the exhaust plume pressure contours on the body and in the sym- 
metry plane. The plume shock is readily evident, as is the footprint of the exhaust flow on the afterbody. 
Because the exhaust flow is at a relatively low Mach number ( M e fa 1 .2 ), expansion waves propagate 
into the plume at a large angle to the local flow direction, and expansion occurs a relatively short distance 
downstream of the nozzle. An expansion fan is seen behind the plume shock, a result of the nozzle flow ex- 
panding along the upswept afterbody. Near the end of the body, a recompression brings the overexpanded 
flow back up to ambient pressure. 

Figure 30 presents surface pressure coefficient comparisons on the windward symmetry plane down- 
stream of the nozzle exit for a series of exit pressure ratios. The data show fairly good correspondence for 
all values of p e /poo- Expansion to ambient pressure takes place within about 2.5 cm (« 5 nozzle heights). 
The calculated pressure is seen consistently to overpredict the data in the initial expansion region of the 
exhaust plume. This is due, at least in part, to nonuniform experimental exit plane conditions wherein the 
exit Mach number (and hence, the pressure) varied circumferentially across the nozzle. Exit plane mea- 
surements show that the exit plane pressure was highest at the symmetry plane. Because the calculations 
assumed uniform exit plane conditions, the flow did not expand as rapidly at the symmetry plane as in the 
experiment. From these comparisons, it is expected that the calculated normal force on the afterbody will 
be slightly higher than what would have been obtained experimentally; a corresponding increment to lift 
and thrust will result. 

With that in mind, lift and thrust are plotted as a function of exit pressure ratio in figure 31. Because 
this model has no inlet, the forces at the nozzle face were adjusted by subtracting the freestream momentum 
flux and pressure. This renders a somewhat realistic representation of the actual forces on the entire vehicle, 
by assuming that flow conditions at the inlet are freestream. The variation of lift and thrust with exit 
pressure is linear, commensurate with the results of section 3.6.1 for the nozzle/afterbody model. The figure 
indicates that net thrust is first generated at p e /Poo & 2.25. Table 3 lists the incremental contributions to 
lift and thrust from the nozzle and afterbody at p e /poo = 6 .27. Incremental values represent the difference 
between the power-on quantity and the power-off quantity. The table indicates that the thrust contribution 
from the afterbody is negligible compared to the contribution from the nozzle face. However, a significant 
lift increment is generated by the afterbody. 
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Table 3. Incremental lift and thrust due to nozzle and afterbody. 
Moo =146, a = 0°, p e /poo = 6.27. 


Component 

cr - ct 0 

o 

0 

1 

Nozzle 

0.1154 

0. 

Afterbody 

0.0049 

0.0151 

Vehicle 

0.1206 

0.0160 


Nozzle and afterbody forces are also important to vehicle balance. As forces on the nozzle face and 
afterbody increase, vehicle balance changes. The pitch-axis moment coefficient, Cm y , provides a convenient 
way to quantify vehicle balance. Since the vehicle weight and center of gravity were not provided, moments 
were calculated about the nose of the body. Figure 32 shows the variation of c mj with exit pressure. The 
pitch moment decreases with increasing pressure, indicating that the nose-down moment due to additional 
afterbody lift overcompensates for the nose-up moment due to additional thrust forces at the exit face as 
the exit pressure increases. Because afterbody lift acts at a large distance from the center of gravity, its 
effect on overall vehicle balance is considerable. Including the moment due to vehicle weight would, of 
course, increase the nominal value of Cm t substantially. 

Limiting streamlines on the afterbody are presented in figure 33 at p e /poo = 6.27. Streamwise flow 
separation is evident near the symmetry plane. Flow leaving the nozzle quickly becomes overexpanded, 
and as it recompresses, the boundary layer separates from the body. The near-surface flow pattern is 
repeated for the lower exit pressure ratio cases, and the separation point moves progressively closer to the 
nozzle as the exit pressure ratio is reduced. This result is only a suggestion of the importance that off- 
design performance analysis may assume. To take a step closer to design-point conditions, though, flow 
solutions are next investigated at a hypersonic Mach number. 

4.6 Generic Hypersonic Vehicle Flow Solutions at Moo - 6 


With a sizeable body of numerical/experimental comparisons in already in hand, a purely numerical 
assay of hypersonic plume/afterbody flow fields is now presented. This set of calculations seeks to: 1) 
demonstrate the ability of the CFD code to solve power-on flow fields in the hypersonic regime, 2) explore 
the fluid dynamics of plume impingement at hypersonic speeds, and 3) investigate the benefits, if any, of 
thrust vectoring for propulsive lift. Regarding the latter, it has been proposed that thrust vectoring may 
be useful at hypersonic flight conditions to trim out very large pitching moments that are generated by 
propulsive forces (private communication from G. S. Deiwert, NASA Ames Research Center, Moffett 
Field, California). 

Selecting flight conditions proved to be somewhat arbitrary. No experimental or flight test data are 
available for powered hypersonic flight as yet. For these numerical experiments, Moo = 6 , J Re = 230 , 000 
cm -1 were chosen to mimic a possible cruise speed (« 17 , 704 m/s) and altitude (« 15 km) of a hypersonic 
transport. Again, fully turbulent flow was assumed over the entire body. Nozzle exit conditions in this 
regime are equally vague, though exit pressure ratios on the order of 100 to 1000 have been spoken of. As 
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such, the operating conditions for these calculations were chosen to be p e /poo = 100, M e = 2 .04. The 
quantity e„ is now introduced to specify the angle of thrust vectoring, measured positive downward from 
the body axis. 

4.6.1 Analysis of the exhaust plume/afterbody interaction- Pressure contours for the entire flow 
field at Moo = 6 are presented in figure 34. The angle of attack, a, is zero, and the thrust vector angle, e„, is 
also zero in this solution. Formation of the wing shock is evident at the wing-body junction, and the plume 
shock is visible on the symmetry plane. The plume shock is clearly defined in figure 35, and its shape 
approximates a conical surface. It interacts very strongly with the afterbody, and figure 36 demonstrates 
the effect it has on flow near the surface. The limiting streamlines on the afterbody reveal a complicated 
flow pattern. Postulated skin friction lines presented alongside the calculated limiting streamlines suggest 
multiple crossflow separations and reattachments, an interpretation of which is put forth here. 

The intersection of the plume shock and afterbody induces crossflow separation S 3 . The thickening 
boundary layer in advance of this separation generates compression waves that coalesce to form a weak 
shock upstream (in the crossflow sense) of the plume shock. The result is a lambda shock structure at the 
body, as depicted in figure 37. Referring back to figure 35, the upstream leg of the lambda shock is evident 
near the body. The upstream leg of the shock induces flow separation along line Si that rolls up into a 
vortex. Unfortunately, tracing particle motions in the vortex does not produce the tightly braided patterns 
seen in subsonic [54] and transonic [55] flows because the streamwise flow speed is so much larger than 
the crossflow speed in hypersonic applications. Nonetheless, the core is evident as a pressure trough in 
figure 37, and a vorticity peak in figure 38. The postulated crossflow pattern described here is sketched in 
figure 39. 

The vortex produces a comer-like flow, driving secondary and tertiary vorticity peaks outboard of it. 
The secondary and tertiary vortices are not actually recirculating, but their influence on the body produces 
the reattachment R\ and separation Si . The large vortex attaches at Ri , just outboard of the strong leg 
of the plume shock. That shock causes flow separation at & , as mentioned before. As the underbody 
fairing tapers, expansion waves radiate from it and bring about reattachment at R 3 . Finally, the symmetry 
plane becomes a separation point by necessity, though the limiting streamlines in figure 36 show only 
slight convergence toward this line. This interpretation of the flow field suggests that the conical afterbody 
fairing is only responsible for R 3 : the other separation and reattachment lines are a direct result of the 
plume shock interaction with the afterbody. Therefore, analogous patterns would be expected for other 
afterbody shapes under similar conditions. 

4.6.2 Effect of thrust vectoring- It has already been observed that most of the scramjet thrust results 
from forces exerted at the exit plane of the nozzle. Because the nozzle is below the body centerline (the 
approximate location of the center of gravity), the thrust force contributes a large nose-up pitching moment. 
Trimming the vehicle with aerodynamic control surfaces becomes impractical; this would require large, 
heavy surfaces and would incur an excessive trim drag penalty. However, the pitching moment due to 
propulsion can be reduced by vectoring thrust from the nozzle downward. As figure 40 shows, thrust 
vectoring rotates the thrust vector toward the center of gravity, effectively reducing the lever arm from 
which this force acts. Thrust vectoring will also produce several potentially important auxiliary effects. 
Lift will increase due to the downward momentum flux from the nozzle, while thrust will decrease. The 
plume shock will be stronger because the external flow must turn more to accomodate the plume flow. 
Finally, the plume flow will expand more as it turns up the afterbody, giving rise to reduced afterbody 
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forces. To determine the net result of these differences, solutions were obtained at e v = 3° and e v = 10°, 
other conditions constant. 

The variation of c my with thrust vector angle is shown in figure 41. As expected, thrust vectoring 
contributes a nose-down pitching moment. In contrast to the sensitivity of c my to changes in exit pressure 
ratio discussed above, the afterbody forces have little effect on moments induced by thrust vectoring. 
Being able to control the vehicle balance in this way may be important to ensure good handling qualities. 
Auxiliary effects of thrust vectoring on lift and thrust are shown in Table 4. At the relatively small vectoring 
angles, the thrust is not seriously degraded, while lift nearly doubles in the interval from e„ = 0° to 
e„ = 10°. 


Table 4. Vehicle lift and thrust with thrust vectoring. 

Moo = 6, a = 0°,p e /poo = 100. 



0° 

3° 

10° 

C T 

0.1173 1 

0.1171 

0.1139 

CL 

0.0156 

0.0202 

0.0307 


Thrust vectoring also deflects flow away from the afterbody. Figure 42 compares limiting streamlines 
on the afterbody for c„ = 0° and e v = 10°. The afterbody flow pattern, though topologically similar to the 
€„ = 0° solution, exhibits less turning of the flow near the body, weakening the exhaust plume/afterbody 
interaction. 


4.7 Summary 


At this point, numerous flow solutions have been presented for different geometries and flow condi- 
tions. Comparisons have been made with experimental data, sensitivities to various parameters have been 
assessed, and interpretations of the flow fields have been made. Any one of these flow solutions contains 
a wealth of information - more than can be analyzed even within the scope of a single paper. The goal 
in this work has been to extract only those aspects of the flow field phenomena that pertain directly to the 
problem at hand. This chapter has addressed those issues in the context of practical flight conditons about a 
realistic geometry. Now it is appropriate to put this work in perspective; to draw more general conclusions 
about the work and its results, and to look for new directions for future work. 
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(a) top view 



(b) side view 


Figure 25. Generic hypersonic vehicle geometry. 
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(c) front three-quarter view 



(d) rear three-quarter view 


Figure 25. Generic hypersonic vehicle geometry, concluded. 
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Figure 26. O-H grid topology for generic hypersonic vehicle. 
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(a) overall view 



(b) detailed view near nozzle exit 

Figure 27. Zonal interface at nozzle exit plane of generic hypersonic vehicle. 
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Figure 28. Perspective view of pressure contours on surface and selected grid 
planes of the generic hypersonic vehicle. Moo — 1.46, a = 0°, 
Pe/Poo = 6.27. 



Figure 29. Detailed view of pressure contours on generic hypersonic vehicle near 
nozzle exit. Moo = 1.46, a = 0°, p e /Poo = 6.27. 
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(b) Pe/Poo = 1.92 



(c) Pe/Poo — 3.97 (d) PtjPoo — 6.27 

Figure 30. Comparison of calculated and experimental surface pressure coeffi- 
cients on windward symmetry plane of generic hypersonic vehicle 
downstream of nozzle exit. M = 1.46, a = 0°. 
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Figure 32. Variation of pitching moment with exit pressure ratio for generic hy- 
personic vehicle. Moo = 1-46, a = 0°. 
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(a) calculated limiting streamlines 



(b) postulated skin friction lines 

Figure 33. Surface flow pattern on generic hypersonic vehicle afterbody. 
Moo = 1.46, a = 0°, Pe/ Poo = 6.27. 
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Figure 34. Perspective view of pressure contours on surface and selected grid 
planes of the generic hypersonic vehicle. M a 0 = 6, a = 0°, 
Pe/Poo — 100. 



Figure 35. Shaded-surface representation of exhaust plume shock for generic hy- 
personic vehicle. = 6, a = 0°, Pe/Poo — 100. 
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(b) postulated skin friction lines 

Figure 36. Surface flow pattern on generic hypersonic vehicle afterbody. 

Moo = 6,a = 0°, pe/Poo = 100. 
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Figure 37. Crossflow pressure contours for generic hypersonic vehicle at 
x = 45 cm. Moo = 6, a = 0°, Pel Poo = 100. 




Figure 38. Crossflow vorticity (z-component of vorticity vector) for generic hy- 
personic vehicle at x = 45 cm. Moo = 6, a = 0°, p e /Poo = 100. 
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Figure 39. Postulated crossflow pattern for generic hypersonic vehicle at 
x — 45 cm. 


81 


'<*> ” '? 3^ 

ay - * T i 1 ». *3 




Bi,. 


iOT FILMED 


EAG£ y ^_.l NI£NIlQNALLX SLM£ 




THRUST 


e 


v 


0 ° 


(a) without thrust vectoring 



(b) with thrust vectoring 

Figure 40. Schematic demonstrating the effect of thrust vectoring on vehicle 
balance. 
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Figure 41. Variation of pitching moment with thrust vector angle 
hypersonic vehicle. M = 6, a = 0°, Pe/Poo = 100. 
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Figure 42. Comparison of calculated limiting streamlines colored by pressure on 
generic hypersonic vehicle afterbody with and without thrust vector- 
ing. M 0 o=6, a = 0°, Pe/Poo = 100. 
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CHAPTER FIVE 


CONCLUSION 


5.1 Summary 


The objective of this investigation was to demonstrate a CFD capability to solve hypersonic flows 
involving exhaust plume/afterbody interaction. The UWIN code was calibrated with all available ex- 
perimental data, and was then extended to conditions of practical interest. Some important parameters 
associated with scramjet propulsion were identified and studied based on the computational results. The 
sensitivity of performance to parametric variations was explained in terms of the underlying fluid dynamic 
phenomena. Some of the first power-on flow solutions have been obtained at hypersonic speeds for a re- 
alistic geometry. The impact of afterbody forces on vehicle balance was assessed, and the complex plume 
expansion process has been studied from a fluid dynamic perspective. 

The nozzle/afterbody model brought out many features characteristic of exhaust plume flows. Com- 
plex flow patterns arise from interference of adjacent axisymmetric nozzles, leading to a highly nonuniform 
pressure distribution and recirculating flow in the footprint region. Variations in afterbody forces with exit 
pressure ratio, upsweep angle, exhaust gas isentropic exponent, and nozzle geometry were assessed, and 
comparisons with experimental data were generally good for these solutions. The experimental compar- 
isons made in this study instilled confidence in the physical correctness of the calculations. Attention 
focused on what more could be learned about the problem from the detailed and voluminous information 
available from a CFD solution. 

Solutions for the generic hypersonic vehicle demonstrated some of the first tip-to-tail flow field solu- 
tions at supersonic and hypersonic Mach numbers for power-on flows. Numerical/experimental compar- 
isons again showed very good correlation. Row field analysis of the supersonic results went on to reveal 
stream wise separation on the afterbody at all exit pressure ratios tested. Afterbody forces at the highest exit 
pressure ratio tested, p e /p<» = 6 .27 , were found to augment gross thrust by 5%, and accounted for 95% 
of total lift. The aerodynamic center was found to move forward 0.5 body lengths from power-off to the 
maximum exit pressure ratio. Predicting such quantities accurately will be critical to designing a vehicle 
with minimal trim drag. In light of the current level of research activity in this field, these results portend 
the significant contribution CFD is now prepared to make to NASP and related hypersonics research. 

Hypersonic results for the generic hypersonic vehicle proved fascinating from a fluid dynamics stand- 
point. The flow conditions prescribed here (Afoo = 6,p e /p <x = 100, Af« = 2. 04), albeit arbitrarily chosen, 
are not unlike those spoken of for a hypersonic passenger transport It was found that the exhaust plume 
shock induces crossflow separation, forming a lambda shock at the body. Vortical flows were found above 
the crossflow separation lines. Within the footprint region, very complicated flow patterns arise from the 
expanding flow and tapering body. Thrust vectoring was shown to weaken the plume/afterbody interaction. 
Because of the magnitude of the thrust, thrust vectoring could be critical to vehicle balance. 
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Finally, the code was shown to be computationally efficient, both for generating solutions from 
freestream initial conditions, and reconverging perturbed solutions (for example, by changing the noz- 
zle exit pressure ratio). A grid resolution study was performed for the nozzle/afterbody model, the results 
of which indicated that the grid dimensions chosen for these calculations were sufficient to ascertain in- 
tegrated quantities, such as afterbody forces and moments, as well as resolve flow field details, such as 
crossflow separation patterns. 

In conclusion, an important first step has been taken in solving the scramjet propulsion problem. 
Inevitably, an introductory investigation such as this raises more questions than it answers, though. Expe- 
rience gained during the course of this investigation indicates several avenues along which further fruitful 
research may be pursued, and these are discussed below. 

5.2 Recommendations for Future Work 


One aspect of interest that is already under investigation concerns the effect of the inlet swallowing of 
some of the external flow on the afterbody flow field (private communication from J. Flores, NASA Ames 
Research Center, Moffett Field, California). In all the cases presented here, the inlet was not modeled. 
This effectively displaces more mass around the body, and one would expect this to result in a higher 
ambient pressure at the nozzle exit. Modeling the inlet with appropriate mass flow will effect a lower 
pressure environment at the nozzle exit, causing the plume to expand more rapidly into the external flow. 
The magnitude of this influence on afterbody forces needs to be determined. 

Whenever discussions of flow separation and skin friction patterns arise, turbulence modeling be- 
comes a central, controversial topic. Kaynak [48] implemented a relaxation modification to the Baldwin- 
Lomax model, the same turbulence model included in the UWIN code, that resulted in much better com- 
parisons between predicted and experimental separation patterns. However, the method addressed only 
transonic shock-induced separation. It may be totally inappropriate for hypersonic flow, and regardless, 
an all-purpose hypersonic turbulence model (ideally including a transition model) would be of far greater 
value. 

Stepwise advancements toward progressively more realistic flight conditions are certain to be forth- 
coming. Enhancements to look for include modeling equilibrium air chemistry for the external flow and 
equilibrium hydrogen-air chemistry for the plume flow. Allowing the streams to mix will also require a 
turbulent shear layer model and a model for turbulent diffusion. Solution-adaptive grids will doubtlessly 
be a necessity for such calculations. Finally, nonequilibrium chemistry in the plume will be necessary to 
predict very high speed, high altitude flight conditions. 

Clearly, CFD research in hypersonic propulsion is still in its infancy. Analysis tools now being devel- 
oped offer fertile ground for technical innovation. And for the first time since the inception of the science, 
computational fluid dynamics stands to play the lead role in the preliminary design of an entire vehicle. 
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APPENDIX 


SIMPLIFYING ASSUMPTIONS FOR THE 
SPECIES CONTINUITY EQUATION 


The full species continuity equation for a binary mixture was presented in equation (2.9). This equa- 
tion was then simplified by neglecting the terms related to diffusion and production. The following pages 
present an order of magnitude analysis of the species continuity equation with application to the present 
problem. This analysis will show that, for the nozzle/afterbody problem addressed here, the diffusion terms 
may be neglected with negligible error, thereby permitting solution of a considerably simpler equation. 

To begin, consider the full species continuity equation for a binary mixture, given in equation (2.9): 

convection <fiffasion 

d(pCt) t d( + ji, ) | d(pcjV + ;,;) + djpcjw + ;,,) 

dt dx dy dz (Al) 

= 

net production 

To arrive at the simplified equation, it must be shown that: 1) there is neither production nor destruction 
of species, and 2) the convection terms are much larger than the diffusion terms. The first condition is 
easily satisfied by noting that SF6 is inert in air, so no chemical reactions take place. Therefore, species 
can neither be produced nor destroyed and K, = 0 . 

To satisfy the second condition, inspect the diffusive mass flux vector, given in equation (2.10) as 

J« = ~~~ rnjmjDjjdj — Df — — i= 1,2,3,- ■■v (A 2) 

p ;-i ° 

where d ; -, the macroscopic gradient vector, is given by 


d = 3(Ty/n) (n L _ 

; dr \ n 


concentraiion 


P 

pressure 


31n(p) 

dr 


-(v-) 




*-1 


force 


(A3) 


The type of diffusion associated with each term is indicated beneath that term. Of these, concentration 
diffusion is the only one that may be important in equation (A.l). The second term, pressure diffusion, 
arises because heavier molecules move more slowly than light ones at the same pressure. In the presence 
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of a pressure gradient, the net diffusion rate of light molecules toward the low pressure is greater than 
that of heavy molecules. Howe and Sheaffer [44] showed that this type of diffusion is important only in 
regions where a pressure gradient coincides with a nonuniform mixture (0 < c < 1 ). Referring to figure 
7, it is seen that for an underexpanded jet flow field, such conditions exist only very close to the nozzle 
lip, where the exhaust gas first expands to ambient pressure. Elsewhere, the mixing layer occurs at the slip 
surface between the plume flow and the external flow. Since the slip surface is by necessity one of constant 
pressure, pressure diffusion is negligible away from the nozzle lip. Force diffusion arises from body forces 
acting nonuniformly on the fluid. In keeping with the assumptions made for the Navier-Stokes equations, 
body forces are neglected and hence the third term in equation (A. 3) vanishes. 


Return now to equation (A.2), inserting the only surviving term from equation (A.3): 

„ 2 2 

dr 


j _ n “ V" m m D 3( n j / K) 
J i / ^ ' 


-D 


p P 


t 


t dln(T) 

j3r 

thermal 


concentration 


(A A) 


Note that v * 2 for the binary mixture assumption applied here. The thermal diffusion term arises from a 
similar mechanism to pressure diffusion: heavy molecules move more slowly than light molecules at the 
same temperature. Hence, in the presence of a temperature gradient, light molecules will migrate more 
rapidly to the lower temperature, effecting thermal diffusion. The nozzle/afterbody experiments employed 
cold jets, wherein the jet stagnation temperature was nearly equal to the freestream stagnation temperature. 
Therefore, temperature gradients in the flow field are small, and thermal diffusion is negligible. 


Now expand die summation of equation (A.4) for the first species of the binary mixture: 


Ji * 


n 2 miDn 


d(n 2 / 

b— 


In) 


(A.5) 


_ r?m\m 2 D\ 2 d(n\/n) 

p dr 

where d( ni /n ) /dr = — d(nz/n)/dr, and Du = 0. Using m = pi /mi, n = p/m, and ci = pi/p, 
equation (A.5) simplifies to 

dc i 


ji = -nm 2 Dn 

r 


Now let 7i77i2 ^ p. giving 


ji as — pD\ 2 


dc i 
dr 


An approximate expression for the binary diffusion coefficient, D \ 2 , is given in reference [42] as 


Dn « 2 .628 • 10 


_ 3V /fy^ 


pa* 


(A. 6) 

(A.l) 

(A8) 


where 


D = binary diffusion coefficient, cm 2 /sec 
T = temperature in Kelvin 
m = molecular weight 


p = pressure in atmospheres 
a = molecular diameter in A 
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This expression is valid only for self-diffusion, so to account for different species, an average mass is used: 
m & mi m 2 / j( mi + m 2 ) . To obtain an order of magnitude estimate for the binary diffusion coefficient, 
the following values are estimated for the nozzle/afterbody experiment: 

T « Too * 100 K 

m 46.18 g/gmole 

P Poo = 0158 atm 

a « 10 A 


. Thus, D \2 « 2 .5 cm 2 /sec = 0(1). 


The concentration gradient is approximated as 


dc 

dr 



(A. 9) 


where 6 is the shear layer thickness. Assuming c\ changes from between 0 to 1 over the thickness of the 
shear layer, this term is of order 6 _1 . Shear layers are jet-like flows, for which the thickness scales with the 
distance (lecture notes from B. J. Cantwell, Stanford University, Stanford, California). Suppose 6 x/de, 
where d t is the nozzle exit diameter. Away from the nozzle exit, then, the shear layer thickness will be on 
the order of 1 to 10 cm within the solution domain, and hence the concentration gradient will be, at most, 
1 cm -1 . Combining this with the approximate binary diffusion coefficient, the diffusive mass flux vector 
is 


ji 


= Di2^-=0(l)0(l) = 0(l) 


(A. 10) 


Referring back to equation (A.1), j\ % appears in the same differential with uc. We already have 
c = 0(1), and if u = O(uoo) = 100,000 cm/sec, then the convective term is 0(100,000). Clearly 
diffusion is negligible where these assumptions hold. 

Further analysis is necessary to determine the role diffusion in the nozzle exit region plays on the 
afterbody forces futher downstream. Another important source of diffusion, that resulting from turbulent 
convection, has been ignored altogether. Finally, in the boundary layer, where the flow velocity vanishes 
at the wall, diffusion also takes on greater importance. But because the pressure is nearly constant in the 
boundary layer, the surface forces are determined by the flow outside the boundary layer, where the fore- 
. going assumptions are valid. Therefore, the diffusion terms are safely omitted for the flow fields presented 
here. 
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